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ABSTRACT 


This  thesis  contains  a  software  package  and  user  documentation  for  use  in 
solving  nonparametric  statistics  problems.  The  following  11  tests  are  considered: 
Binomial  test,  Sign  test,  Cox-  Stuart  test,  Quantile  test,  Quantile  Confidence  Intervals, 
McNemar  test  (N  £  20),  Tolerance  Limits,  Cochran  test,  Chi  Square  test,  Median 
Test,  and  Chi  Square  Goodness  of  fit  tests.  The  programs  are  in  BASIC  and  are 
designed  specifically  for  use  on  the  Tandy  Radio  Shack  M 100  microcomputer. 


Accesion  For 

J 

NTIS  CRA&I 

« 

DTIC  TAB 

o 

Unannounced 

u 

Justification . . 

- 

By  .  .. 

Disti  ibution  / 

Availability  Codes 

Dist 


Avail  and  I  or 
Special 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 

II.  PROGRAM  DESIGN  ISSUES . 

A.  EQUIPMENT  AND  LANGUAGE  . 

B.  KNOWLEDGE  LEVEL  OF  THE  USER . 

C.  SELECTION  OF  TESTS . 

D.  SCREEN  DISPLAYS . 

E.  DESIGN  OF  THE  USER'S  MANUAL . 

III.  DISCUSSION  OF  THE  PROGRAM  MODULES . 

A.  GENERAL . 

B.  PROGRAMS  . 

1.  Binomial  Modules  . 

2.  Chi  Square  Module . 

IV.  AREAS  FOR  FURTHER  WORK . 

APPENDIX  A:  DOWNLOADING  INSTRUCTIONS . 

APPENDIX  B:  USER  S  MANUAL . 

1.  GENERAL  INFORMATION . 

2.  TESTS  BASED  ON  THE  BINOMIAL  DISTRIBUTION  . 

a.  Binomial  Test . 

b.  Quantile  Test . 

c.  Sign  Test . 

d.  Cox-Stuart  Test . 

e.  Quantile  Confidence  Interval  . 

3.  TESTS  BASED  ON  THE  CHI  SQUARE  DISTRIBUTION 

a.  Tolerance  Limits  . 

b.  Cochran  Test . 

c.  Median  Test . 


j 


( 


d.  Chi  Square  Test  for  Differences  in  Probabilities 

e.  McNemar  Test . 

f.  Goodness  of  Fit  test . 

APPENDIX  C:  MODIFYING  PROGRAM  MODULES . 

APPENDIX  D:  POSSIBLE  PROGRAMMING  PROBLEMS  . 
APPENDIX  E:  BINOMIAL  DISTRIBUTION  PROGRAM  . . 
APPENDIX  F:  CHI  SQUARE  DISTRIBUTION  PROGRAM 

LIST  OF  REFERENCES . 

BIBLIOGRAPHY . 

INITIAL  DISTRIBUTION  LIST . 


I.  INTRODUCTION 


Nonparametric  statistics  is  an  area  in  which  mastery  requires  the  solution  of 
numerous  time  consuming  problem  sets.  Because  of  this,  students  sometimes  do  not 
master  the  techniques.  Thus  the  first  purpose  of  this  thesis  is  to  provide  a  software 
package  to  assist  students  in  solving  problems  involving  eleven  nonparametric  tests. 

The  programs  are  written  for  the  Tandy  Radio  Shack  M100,  a  laptop  computer 
with  24k  RAM.  These  computers,  complete  with  a  communications  cable,  an  AC 
adapter,  and  an  owner's  manual  are  issued  to  students  in  the  Operations  Research 
Department  at  the  Naval  Postgraduate  School.  However,  for  various  reasons,  they  are 
underutilized.  This  thesis  will  provide  the  students  with  programs  which  are  directly 
applicable  to  course  work,  thus  improving  utilization. 
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II.  PROGRAM  DESIGN  ISSUES 

This  Chapter  presents  a  brief  overview  of  the  design  decisions  made  in  developing 
the  software. 

A.  EQUIPMENT  AND  LANGUAGE 

The  M100  computer  was  selected,  since  this  is  the  laptop  portable  issued  to  the 
students.  The  software  package  was  written  in  BASIC  for  several  reasons.  First, 
BASIC  is  the  ROM  resident  language  available  on  the  M100.  This  was  the  most 
important  reason  for  its  selection.  Also,  the  most  common  microcomputer  language  is 
BASIC  [Ref.  l:p.  185].  Although  the  modules  were  written  in  a  BASIC  specific  to  the 
M100,  they  could  be  modified  with  relative  ease  to  run  on  other  microcomputers.  The 
use  of  assembly  language  would  not  have  made  the  programs  significantly  smaller, 
since  over  50%  of  the  program  is  used  for  screen  displays. 

B.  KNOWLEDGE  LEVEL  OF  THE  USER 

It  was  assumed  that  the  primary  user  of  the  software  package  will  have  had  the 
three  course  series  in  parametric  statistics  taught  at  the  Naval  Postgraduate  School, 
and  will  either  have  completed,  or  be  currently  enrolled  in  a  nonparametrics  course.  It 
is  also  assumed  that  the  user  will  have  had  exposure  to  the  IBM  370  and 
microcomputers  before  attempting  to  use  the  programs. 

C.  SELECTION  OF  TESTS 

This  software  package  was  designed  specifically  to  accompany  and  parallel  the 
main  statistical  tests  covered  in  chapters  three  and  four  of  the  Conover  text  on 
nonparametrics  used  at  the  Naval  Postgraduate  School  [Ref.  2:pp.  95  -  212].  It  is 
hoped  that  the  programs  will  be  immediately  applicable  to  course  work  being  taught 
here.  Not  all  tests  from  these  two  chapters  were  included,  because  the  limited  memory 
capacity  of  the  M100  restricts  the  number  of  programs  which  could  be  included  and 
still  permit  an  easy  transfer  of  programs  from  the  mainframe  via  the  modem. 
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D.  SCREEN  DISPLAYS 


Some  of  the  design  concepts  used  in  programs  for  Computer  Assisted  Instruction 
were  employed  in  the  screen  displays  and  module  structures.  The  interactive  screen 
displays  were  designed  to  avoid  excessive  verbiage  and  an  attempt  was  made  to  vary 
screen  displays  in  order  to  maintain  the  interest  of  the  user  [Ref.  3:p.  29].  To  provide 
the  user  with  the  feeling  of  control  over  the  programs,  all  the  screen  displays  remain 
available  for  viewing  until  input  is  received  from  the  user  to  indicate  that  he  wishes  to 
proceed. 

The  modules  were  designed  as  single  level  menu  driven  programs.  This  is 
appropriate  because  the  software  package  is  designed  to  accompany  a  textbook  as  a 
learning  aid  for  students.  In  addition,  multi-level  menus  require  additional  RAM 
which  is  not  available  on  the  OR  department's  M100  machines. 

The  problem  solving  approaches  used  in  the  programs  follow  those  in  the 
textbook.  The  same  logic  patterns  were  employed  in  the  approach  to  each  problem  in 
order  to  improve  the  problem  solving  techniques  of  the  students  using  the  programs 
[Ref.  l:pp.  298  -  305], 

E.  DESIGN  OF  THE  USER'S  MANUAL 

The  user's  documentation  is  divided  into  four  sections.  The  first  contains  the 
information  needed  by  a  novice  user  to  download  the  programs  from  the  mainframe 
using  the  Ml 00  (See  Appendix  A).  Next  are  the  general  instructions  on  the  module 
programs,  and  example  problems  of  each  test  with  screen  displays  shown.  Most  of  the 
example  problems  included  demonstrate  the  application  of  nonparametric  statistics  to 
military’  problems.  This  portion  of  the  documentation  is  included  as  the  User  s  Manual 
in  Appendix  B.  The  final  two  portions  of  the  documentation  cover  possible  problems 
with  the  programs  and  procedures  for  redimensiomng  arrays.  (See  Appendix  C  and  D). 


III.  DISCUSSION  OF  THE  PROGRAM  MODULES 

A.  GENERAL 

The  package  has  5  subsections.  The  first  contains  information  on  downloading 
programs  from  the  mainframe,  a  user's  manual  for  the  programs,  instructions  on 
redimensioning  the  program  arrays,  and  possible  solutions  to  program  run  errors.  This 
document  is  designated  CONOVER  SCRIPT,  and  is  available  on  the  mainframe 
"M100"  disk.  This  subsection  is  designed  to  be  printed  using  the  Naval  Postgraduate 
School  laser  printer  from  a  3278  terminal. 

The  last  four  components  are  the  actual  programs  themselves.  These  routines 
are  stored  as  BINOMA  SCRIPT,  BINOMB  SCRIPT  for  the  tests  covered  in  chapter  3 
of  Conover  and  CHISQA  SCRIPT  and  CHISQB  SCRIPT  for  the  tests  derived  from 
the  chi  square  distribution  covered  in  chapters  3  and  4  of  Conover.  [Ref.  2:pp.  95  - 
212).  These  modules  are  also  located  on  the  "M100"  disk. 

B.  PROGRAMS 

1 .  Binomial  Modules 

BINOMA  and  BINOMB  compose  the  two  subpackages  which  must  be 
downloaded  in  order  to  calculate  the  solutions  to  problems  using  the  following 
procedures: 

1.  Binomial  Test, 

2.  Quantile  Test, 

3.  Sign  Test, 

4.  Cox-Stuart  Test, 

5.  Quantile  Confidence  Interval. 

The  user  is  prompted  to  input  values  needed  in  order  to  calculate  the  test  statistic,  and 
is  also  permitted  to  select  one-sided  or  two-sided  test.  The  user  may  also  select  any 
type  I  error  rate  (Alpha)  desired  in  conducting  the  hypothesis  test.  The  Alpha  value 
returned  by  the  machine  is  rounded  down  if  the  user's  alpha  cannot  be  determined 
exactly,  due  to  discreteness.  After  the  calculations  of  the  hypothesis  test,  the  user  is 
given  the  opportunity  to  calculate  the  actual  significance  level  (Alpha  Mat)  based  on 
the  sample  data,  return  to  perform  another  test,  or  breakout  to  the  M100  main  menu. 
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2.  Chi  Square  Module 

The  CHISQA  and  CHISQB  cover  six  procedures  which  use  the  chi  square 
distribution.  These  are  covered  in  Chapters  3  and  4  of  the  Conover  text  [Ref.  2:pp.  95  - 
212].  These  are  as  follows: 

1.  Tolerance  Limits, 

2.  McNemar  Test, 

3.  Cochran  Test, 

4.  Median  Test, 

5.  Chi  Square  Test, 

6.  Goodness  of  Fit  Test. 

Unlike  the  Binomial  programs,  these  procedures  require  sample  data  for  the 
performance  of  calculations.  Yates'  correction  for  continuity  is  not  included,  since  this 
is  Conover's  recommendation  [Ref.  2:p.  149].  At  the  conclusion  of  the  calculations  for 
each  test,  the  user  is  provided  the  value  of  the  test  statistic  and  the  alpha-hat  value 
(level  of  significance)  associated  with  that  test  statistic.  The  user  can  determine  if  the 
alpha-hat  value  is  small  enough  for  the  rejection  of  the  null  hypothesis.  In  the  case  of 
Goodness  of  Fit  tests,  the  program  tests  only  normal,  poisson,  and  uniform 
hypotheses.  In  both  the  normal  and  the  poisson  cases,  the  degrees  of  freedom  are 
adjusted  if  parameters  are  estimated  from  the  sample  data. 


IV.  AREAS  FOR  FURTHER  WORK 


There  are  several  possible  extensions  of  this  work.  It  would  be  useful  to  have 
programs  for  the  tests  in  Chapters  5  and  6  of  the  Conover  text,  on  additional  modules 
[Ref.  2:pp.  213  -  385].  This  would  provide  programs  which  would  parallel  the  entire 
course  in  nonparametrics  as  it  is  taught  at  the  Naval  Postgraduate  School. 

If  the  RAM  capacity  were  upgraded  to  32k  on  the  M100  machines,  the  programs 
could  be  collapsed  into  fewer  modules.  This  would  simplify  the  transfer  of  the 
programs  via  modem. 


APPENDIX  A 

DOWNLOADING  INSTRUCTIONS 

General 

The  following  instructions  are  provided  in  addition  to  those  included  in  Tandy 
Owner's  Manual  [Ref.  4:pp.  75  -96],  to  assist  you  in  interfacing  with  the  IBM  370  at 
the  Naval  Postgraduate  School,  and  to  more  efficiently  download  the  "BINOM"  and 
the  "CHISQ"  modules  located  on  the  M100  library.  If  you  have  previously 
downloaded  programs,  and  your  A-disk  contains  the  necessary  files,  you  will  wish  to 
skip  the  first  section  of  these  instructions  and  go  directly  to  the  instructions  on  the 
modules  themselves. 


Preparing  Your  A- Disk 

Before  attempting  to  download  any  programs  from  the  M100  library,  you  should 
perform  the  following  actions  from  a  standard  3278  terminal  [Ref.  5:p.  1]. 

1.  Add  a  LINKMIOO  EXEC  to  your  A-disk.  It  should  contain  the  following  lines: 

CP  LINK  0344p  195  295  M100 
ACC  295  B 

2.  From  the  CMS  mode,  link  to  the  M100  library  using  the  LINKMIOO  command. 
You  should  either  create  a  Profile  Exec  identical  to  the  the  Profile  Exec  located 
in  the  M 100  library,  or  augment  your  present  Profile  Exec  with  the  instructions 
on  the  M100  Profile  Exec  [Ref.  5:p.  I).  The  simplest  way  to  do  this  is  to  enter 
the  Xedit  Profile  Exec  and  use  the  command  GET  PROFILE  EXEC  BI. 

If  you  expect  to  upload  files,  you  should  create  or  augment  your  Profile  Xedit  with  the 
Profile  Xedit  information  from  the  M100  library  in  the  same  fashion  that  you  modified 
your  Profile  Exec.  (Note  -  be  certain  that  your  Profile  contains  Set  Linedit  off  in  order 
to  prevent  data  transfer  problems.) 

While  at  your  terminal,  if  you  would  like  to  examine  any  of  the  files  located  in  the 
M100  library,  simply  type  flist  *  *  B  and  use  the  "Browse''  command,  since  this  is  a 
read  only  disk.  These  are  the  only  changes  to  your  A-disk  which  arc  necessary  to  use 
your  M100  from  home  to  download  M100  library  programs. 
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Downloading  Procedures  Using  the  TELCO M  Program 


Initially,  it  will  be  necessary  to  use  the  ROM  resident  TELCOM  package  on  the 
M100  to  download  the  CMSCOMMS  program.  This  program,  written  by  James  Eagle, 
simplifies  the  M100  interface  with  the  IBM  370.  This  program  will  be  used  to 
download  the  program  modules  from  the  mainframe. 

You  will  find  the  instructions  on  setting  the  switches,  and  connecting  your  M100 
with  the  modem  in  the  Tandy  Owner's  Manual  [Ref.  4:pp.  76-77).  After  connecting  the 
modem  cable  to  the  MLOO  and  your  telephone,  use  the  following  procedures  to 
download  programs  using  TELCOM. 

STEP  1A: 

Turn  on  the  M100  and  access  the  TELCOM  program.  (The  " label"  keys  above  the 
keyboard  will  automatically  display  the  assigned  function  key  labels)  [Ref.  4:p.  79] 

STEP  2A: 

Strike  the  "J2"  key  function  key.  The  word  "CALL"  will  be  displayed.  Enter  the 
following  numbers  and  symbols  exactly  as  they  appear  here: 

6463025<> 

(The  extra  symbols  ( <  >  )  are  required.)  Strike  the  "enter”  key.  The  M 100  will  dial  the 
number.  When  the  connection  is  being  made  with  the  IBM,  you  will  hear  an  audible 
beep.  (If  you  do  not  get  a  beep,  the  M100  you  have  has  had  the  command  "Sound  OfT 
input.  If  you  want  to  turn  it  back  on,  see  the  Owner's  Manual)  [Ref.  4:p.  123].  When 
the  connection  is  established,  the  beeping  stops  and  the  the  screen  displays  "VM;370 
ONLINE'. 

Strike  the  "enter"  key.  An  exclamation  mark  (!)  will  appear  on  the  left  side  of 
your  screen,  and  a  period  (.)  will  appear  on  the  right.  The  period  (.)  is  the  "ready" 
prompt.  You  may  now  logon  the  system  using  normal  procedures.  (Note 
transmissions  from  the  mainframe  are  split  between  lines,  with  the  first  character  of 
each  line  displayed  at  the  right  side  of  the  screen.  The  "CMSCOMMS"  software 
package  will  eliminates  this  inconvenience.) 
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When  the  machine  queries  you  for  your  password,  first  a  row  of"***"  are  printed, 
then  these  symbols  are  overprinted  by  a  series  of  "HHH".  Finally  the  machine  will 
print  a  series  of  "SSS"  proceeded  with  a  period  (.)-  You  should  not  attempt  to  input 
your  password  until  the  "SSS"  series  has  completed  printing.  If  you  receive  an  error 
message,  reinitiate  the  logon  procedure.  (Note  -  the  mainframe  will  transmit  an  error 
when  you  input  a  line  in  which  you  have  backspaced  to  correct  a  character.) 

STEP  3 A: 

The  word  “Half  should  be  displayed  above  th e"/4"  key.  If  it  is  not,  strike  the  ”f4“ 
function  key  once. 

STEP  4 A: 

If  the  screen  has  not.  yet  indicated  that  “CMS"  is  on  line,  enter  either  I  CMS,  or  B. 

STEP  5 A: 

Link  to  the  M100  library  by  typing  LIN  KM  100.  If  you  would  like  to  examine  the 
CMSCOMMS  file  before  downloading  it,  input  TYPE  CMSCOMMS  JNE  BI,  and 
strike  the  "enter"  key.  This  allows  you  to  "Browse"  the  file  as  it  scrolls  on  the  terminal 
screen.  The  "W A  IT'  command  will  be  displayed  on  the  bottom  of  the  screen  above 
the  "f6"  key  while  the  file  scrolls.  When  the  end  of  file  is  reached,  the  “WAIT'  display 
disappears,  and  the  machine  will  return  an  " R "  on  the  right  side  of  the  screen.  This 
indicates  that  you  are  now  in  CMS. 

STEP  6A: 

Strike  the  "f2"  function  key  for  "down".  The  M100  will  query  you  “File  to  be 
downloaded!" .  You  may  type  in  any  6  letters  you  wish  to  use  to  designate  the  file  name 
to  be  used  in  the  M100  RAM  file  directory.  It  is  suggested  that  you  select  the  name 
CMSCOM,  since  we  have  used  this  name  throughout  these  directions.  Type  the  name 
and  strike  the  "enter"  key.  This  will  cause  the  "down"  label  on  the  bottom  of  the  Ml  00 
display  to  change  to  reverse  video.  We  are  now  ready  to  transfer  the  file. 

STEP  7A: 

Type  in  TYPE  CMSCOMMS  JNE  BI,  and  strike  the  "enter"  key.  The  program 
you  saw  in  Step  5  should  scroll  on  the  screen.  If  you  accidentally  typed  the  line  before 
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the  "ready"  prompt  was  displayed,  you  will  receive  an  error.  Wait  for  the  return  of  the 
prompt  (.)  and  reenter  your  line.  If  your  download  is  uneventful,  wait  until  you  receive 
an  "R"  from  the  IBM,  indicating  it  has  returned  to  CMS.  This  "R"  will  be  displayed  on 
the  right  side  of  the  screen.  Strike  the  “JT  for  the  " down  function  key  again.  The 
"down”  which  was  previously  displayed  in  reverse  video,  will  now  go  back  to  normal 
video. 

STEP  8A: 

The  transfer  is  complete  and  you  may  logoff  the  mainframe.  Type  LOGOFF  and 
the  IBM  will  return  the  session  statistics.  When  it  has  completed,  input  "BYE'  by 
striking  the  "f8“  key.  The  machine  will  query  with  "disconnect!".  Respond  with  YES, 
or  Y,  and  strike  the  "enter"  key.  The  function  keys  labels  will  change.  We  return  to  the 
M100  main  "menu",  by  striking  the  “f8“ key. 

Converting  An  ASCII  File  to  a  Basic  Program  on  the  MI00 

The  file  "CMSCOM.DO"  is  now  displayed  on  the  M100  main  menu.  The  ".DO" 
indicates  that  the  file  is  in  ASCII  or  character  format.  To  run  the  program  as  a  BASIC 
program,  we  must  use  the  text  editor  to  eliminate  the  heading  and  trailing  "garbage"  • 
associated  with  the  data  transfer,  and  remove  unnecessary  carriage  returns.  IF  YOU 

ARE  INEXPERIENCED  WITH  THE  TEXT  EDITOR  REFER  TO  THE  TANDY 
OWNER  MANUAL  PAGES  43  -  58.  Pages  53  and  54  deal  specifically  with  deleting 
text. 

STEP  IB: 

Enter  the  CMSCOM  document  file  and  remove  heading  material,  using  the  text 
editor,  until  the  first  numbered  program  line  showing  reads  as  follows: 

1  CLS:  PRINT"  Terminal  Package  for 

CMS": PRINT 

Strike  <  control  key>  and  the  down  arrow  above  the  keyboard  to  go  to  the  end  of  file. 
Here  remove  trailing  material  until  the  last  line  reads: 

702  CLOSE: MAXFILES=1: MENU 
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STEP  2B: 


Leave  the  text  editor  by  striking  "f8".  Now  enter  the  BASIC  program  and  type  the 
following  command  exactly  as  shown: 

LOAD  "CMSCOM. DO"  <enter  key> 

The  word  "Wait"  will  flash  off  and  on  in  reverse  video  as  the  program  is  being 
compiled.  When  the  M100  returns  an  "OK",  it  means  your  program  has  loaded 
correctly.  (If  instead  of  an  “OK’  after  the  flashing  "Wait”,  you  receive  an  error 
message,  this  indicates  that  you  have  some  carriage  returns  (the  carat  symbols)  at 
locations  other  than  at  the  end  of  lines.  You  will  need  to  return  to  the 
"CMSCOM. DO"  and  scroll  through  it  using  the  text  editor  to  eliminate  unnecessary 
carriage  returns.)  "SAVE'  the  BASIC  program  as  a  RAM  file,  by  typing  in  exactly  the 
following: 

SAVE  "CMSCOM"  <enter  key> 

The  M100  will  respond  with  another  "OK".  You  may  now  strike  "JS"  and  return  to  the 
main  menu.  Both  the  document  program  and  the  basic  program  will  be  shown  as 
RAM  programs.  Erase  the  document  program  by  entering  BASIC  and  typing  the 
following: 

KILL  "CMSCOM. DO"  <enter  key> 


You  now  have  the  program  "CMSCOM. BA"  in  your  M100  and  can  proceed  to 
download  the  program  modules. 

Note  that  you  may  wish  to  also  download  the  program  CROUT  JNE  B1  from  the 
M100  library.  This  can  be  used  to  remove  "most"  of  the  extra  carriage  returns  from 
downloaded  ASCII  files.  You  will  not  need  to  use  this  program  with  the  program 
modules,  because  they  have  been  designed  with  lines  of  80  characters  or  less. 


Downloading  the  BINOM  and  CIIISQ  Program  Modules 
Using  the  CMSCOM.BA  Program 

BEFORE  ATTEMPTING  TO  DOWNLOAD  MODULES  ERASE  ALL  PROGRAMS 
EXCEPT  CMSCOM.BA,  and  CROUT.BA.  Failure  to  do  so  will  cause  your  machine  to 
exceed  US'  24k  memory  capability. 

Since  interfacing  with  the  N’PS  mainframe  is  simplified  by  using  the  CMSCOM.BA 
program  written  by  James  Eagle,  we  will  use  it  to  download  the  program  modules. 
The  steps  to  use  are  as  follows: 

STEP  1C. 

Run  the  " CMSCOM.BA "  program  by  placing  the  cursor  over  this  prog'am.  On  the 
first  screen,  the  program  will  ask  “On-campus  call(y  n)7'.  Answer  .V,  unless  \ou  using  a 
military  phone  on  base.  The  M100  will  dial  the  number.  You  will  hear  a  beep  or 
"ring"  as  long  as  the  system  is  trying  to  establish  a  connection.  When  the  connection 
has  been  made  with  the  mainframe,  the  "beeping"  will  cease,  and  you  should  sa  ke  the 
"enter"  key. 

STEP  2C: 

The  following  data  will  be  added  to  the  screen. 

1463025 

The  period  (.),  as  with  TELCOM,  is  the  "ready"  prompt. 

The  machine  is  now  ready  for  the  logon  procedures.  These  are  identical  to  those 
used  when  you  are  at  the  standard  3278  terminal.  After  you  have  logged  on,  strike  the 
"enter"  key  again  and  wait  for  the  machine  to  return  an  R  before  continuing.  Strike 
the  "enter"  key  once  more.  This  may  take  some  time  if  the  load  on  the  mainframe  is 
heavy.  If  you  had  a  full  screen  shown,  you  woulJ  sec  the  Rl'.XX I. XG  down  in  the 
lower  right  hand  corner  of  the  terminal  screen. 


STEP  3C: 

Strike  the  " label '  key  above  the  keyboard  to  display  the  function  key  assignments 
used  in  the  program.  You  must  now  type  LINKM100,  to  access  the  M100  library 
where  the  modules  are  stored. 

STEP  4C: 

Strike  the"/2"  key  for  "down"  to  begin  downloading  the  program  module.  The 
screen  will  prompt  you  with  "File  to  download ?".  Enter  the  6  letter  designator  you  wish 
to  use  for  the  M100  file.  It  is  recommended  that  if  you  are  unfamiliar  with  these 
procedures,  that  you  use  the  module  names  assigned  in  these  instructions  to  prevent 
confusion.  The  M 100  screen  will  now  display  the  following  line: 

Downloading  to  <6  letter  t1tle>  until  R  entered. 


The  two  programs,  BINOM  and  CIIISQ  are  divided  into  parts 
the  M100.  This  has  been  done  because  the  files,  in  their  entirety, 
novice  user  of  the  M100  to  convert  into  BASIC  easily.  If  you  are 
distribution,  you  must  transfer  both  the  BINOMA  SCRIPT  BI, 
SCRIPT  Bl  files.  Tor  the  Chi  Square  Distribution  problems, 
SCRIP  I  Bl,  and  CHISQB  SCRIPT  Bl.  We  will  demonstrate  the 
two  BINOM  modules.  (Note  -  You  must  download  both  halves 
convert  them  to  BASIC  in  order  to  be  able  to  use  the  programs.) 

STEP  5C: 


for  data  transfer  to 
are  too  large  for  a 
using  the  binomial 
and  the  BINOMB 
transfer  CIIISQA 
download  with  the 
one  at  a  time  and 


Strike  the  "enter"  key  to  bring  up  the  "ready"  prompt  (.).  Next  type  in  TYPE 
BISOMA  SCRIPT  BI  <  enter  key  > .  The  program  you  arc  downloading  should  begin 
to  scroll  by.  If  there  are  any  problems  and  the  program  does  not  begin,  strike  <  control 
key  >  R,  and  respond  to  the  screen  display  "Press  any  key  for  Terminal  Mode".  After 
this  go  to  STEP  4C  and  try  again.  The  program  modules  arc  relatively  long,  so  they 
will  take  approximately  .5-7  minutes  to  complete  the  data  transmission. 


STEP  6C: 


The  receipt  of  an  " R“ ,  means  we  are  now  in  the  CMS  mode  and  that  we  can  input 
<  control  key>  R,  and  Press  any  key  to  return  to  the  terminal  mode. 

STEP  1C: 

Wait  for  the  "ready"  prompt,  then  LOGOFF.  The  M100  will  return  a  dotted 
cursor.  Strike  the  "f8"  to  disconnect  and  return  to  computer  mode. 

Module  Conversion  to  BASIC 

These  instructions  are  very  similar  to  those  used  to  convert  the  CMSCOM.BA 
program  presented,  except  when  the  second  module  is  converted,  and  "pasted"  onto  the 
first  submodule.  These  instructions  should  be  followed  closely  or  you  may  get  an  OM 
(Out  of  Memory)  error. 

STEP  ID: 

Enter  the  BINOMA  document  file,  and  with  the  text  editor  remove  the 
extemporaneous  heading  and  trailing  information  which  will  prevent  its  operation  as  a 
BASIC  program.  When  you  complete  this  operation,  the  first  and  last  lines 
respectively,  should  read  as  follows: 

100  CLS 

14000  '****  CALCULATION  OF  ACTUAL  ALPHA  HATS 


Leave  the  text  editor  and  enter  the  BASIC  program  by  striking  "fS"  for  the  "menu" 
and  then  the  <  enter  key>  .  In  BASIC  type  the  following  command. 

LOAD  "BIN0MA. DO"  <enter  key> 

The  word  "WAIT"  will  flash  off  and  on  in  reverse  video  as  the  program  compiles. 
When  the  M100  returns  an  "OK",  it  means  your  program  has  loaded  correctly.  (An 
error  message,  indicates  that  you  have  some  unnecessary  carriage  returns  in  the 
document.  You  will  need  to  return  to  the  "BINOM.DO"  and  scroll  through  it  using 
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the  text  editor  to  eliminate  unnecessary  carriage  returns.  This  should  not  happen 
unless  you  had  transmission  errors.)  "SAVE"  the  BASIC  program  as  a  RAM  file,  by 
typing  in  the  following: 

SAVE  "BINOMA"  <enter  key> 

The  M100  will  respond  with  another  "OK".  Strike  “J8“  and  return  to  the  main  menu. 
Check  to  ensure  that  you  have  both  the  document  program  and  the  basic  program 
shown  as  RAM  programs.  Enter  the  BINOMA. BA.  The  program  should  immediately 
begin  to  run.  To  escape  from  the  program  hold  down  the  <  shift  key>  and  strike  the 
<  break  key>  above  the  keyboard,  or  use  <  control  key>  C  ,  followed  by  "f8”  to  leave 
the  program.  If  you  have  an  error,  examine  the  program  using  either  the  "LIST" 
command  [Ref.  4:pp.  152-153],  or  the  text  editor.  If  the  program  executed,  however, 
erase  the  document  version  by  entering  BASIC  and  typing  the  following: 

KILL  "BINOMA. DO"  <enter  key> 

To  download  the  next  submodule  repeat  STEPS  Cl  through  C7.  After  the 
submodule  BINOMB.DO  is  downloaded,  follow  the  procedures  listed  below  to  add  it 
to  the  BINOMA. BA  program. 

STEP  2D: 

Enter  the  BINOMB.DO  file  and  eliminate  the  unnecessary  information.  The  first 
and  last  lines  should  be  as  listed  below 

14020  ON  QT  GOTO  14040,14300,14040,13100,14040 

18770  PRINTTAB(5)  "reject  your  null  hypothesis  if": RETURN 

STEP  3D: 

We  will  now  "paste”  the  two  files  together. 

1.  Move  the  cursor  to  the  top  of  the  BINOM.DO  file,  and 
strike  <  enter  key  >  . 

2.  Turn  the  "labels"  on,  using  the  “label"  button. 


3.  Strike  “J7“  to  "select"  text  for  copying. 

4.  Strike  <  control  key  >  and  hold  while 
striking  the  <  down  >  cursor  arrow. 

The  cursor  will  go  to  the  end  of  the  file,  and  the  entire  file  will  be  highlighted  in  reverse 
video. 

5.  Strike  the  "f6"  to  cut  the  program  to  the 
"paste"  buffer  of  the  M 100. 

6.  Strike  "f8",  wait  for  an  "OK",  then  strike  "f8"  again  to 

return  to  the  main  menu. 

7.  Enter  the  "BINOMA. BA"  program,  then  break  out  by  either 

using  <  control  key>  C  or  <  shift  key>  <  break  key>  . 

An  "OK"  will  appear  on  the  screen. 

8.  Strike  the  “Paste"  key  above  the  keyboard,  and  the  new 

program  will  be  " pasted "  to  the  previous  program. 

The  "paste"  buffer  on  the  M100  will  still  contain  the  entire  submodule.  This  will  use 
an  excessive  amount  of  memory  so  you  will  want  to  clear  the  buffer  by  copying 
something  short  into  it.  (It  remains  full  even  when  the  machine  is  turned  off.)  To  clear 
the  buffer,  use  the  Copy  and  Paste  commands  listed  in  the  Tandy  Owner's  Manual 
[Ref.  4:p.  56]. 

You  will  now  need  to  attempt  to  run  the  program  to  determine  if  the  load  has 
been  successful.  A  " DS"  error,  indicates  that  you  have  an  extra  carriage  return  in  the 
program  which  must  be  removed  before  it  can  operate  correctly.  You  will  need  to  edit 
the  program  to  remove  it. 

This  completes  the  transfer  of  the  Binomial  Distribution  program  module.  To  the 
download  the  Chi  Square  Distribution  program  module,  simply  load  CHISQA  SCRIPT 
first,  and  then  CHISQB  SCRIPT  (the  order  is  important  because  of  memory 
capacities). 


REPRODUCED  AT  GOVERN  MEN  I  F'<rnr>b' 
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APPENDIX  B 


USER'S  MANUAL 


1.  GENERAL  INFORMATION 


It  is  assumed  that  users  of  this  software  package  are  familiar  with  nonparametric 
statistics  and  computers.  Since  the  programs  are  written  for  the  TRS  M100  and  use 
the  IBM  370  as  a  storage  device,  a  knowledge  of  micros  or  mainframes  with  virtual 
machine  interfaces  is  required  to  use  the  programs.  If  you  have  not  yet  downloaded  a 
package  module,  see  Appendix  B  on  Downloading  Programs  before  proceeding.  The 
programs  themselves  have  been  especially  designed  to  accompany  W.J.  Conover’s  book 
Practical  Nonparametric  Statistics  (Ref.  2}.  The  symbols,  notation,  and  organization 
used  in  the  programs  parallel  this  text  very  closely. 


Before  you  attempt  to  use  any  of  the  programs,  you  should  ensure  that 


The  "Caps  Lock"  key  is  depressed  * 


on  the  M100.  If  you  fail  to  do  this,  the  M100  will  not  respond  properly  to  your 
keyboard  inputs.  You  should  also  be  aware  of  the  specific  meanings  of  phrases  used 
on  the  interactive  screens. 


1.  "Strike  any  Key"  means  keyboard  keys  only,  not  function  keys, 
shift  keys  or  control  keys. 


2.  "MAKE  YOUR  SELECTION"  requires  the  input  of  an  integer  number 
from  the  keyboard,  not  a  function  key,  shift  key,  or  control 
key. 


If  you  strike  the  keys  too  forcefully,  the  machine  may  skip  screens. 


This  package  is  composed  of  two  software  modules,  "BINOM",  and  "CHISQ".  The 
BINOM"  module  covers  tests  from  Chapter  3  of  Conover  which  have  the  binomial  as 
their  underlying  probability  distribution  function.  The  "CH1SQ"  module  includes  tests 
covered  in  Chapters  3  and  4  which  are  derived  from  the  Chi  Square  distribution. 
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2.  TESTS  BASED  ON  THE  BINOMIAL  DISTRIBUTION 


Included  on  the  module  "BINOM"  are  the  binonual  test,  the  quantile  test  and  the 
quantile  confidence  interval,  the  Cox-Stuart  test,  and  the  sign  test.  Since  these  tests 
are  all  derived  from  the  two  parameter  binomial  PDF.  they  all  require  the  input  of  an 
"N",  and  a  "P  "N"  is  equal  to  the  number  of  trials,  and  "P  "  is  the  expected  or 
hypothesized  probability  of  a  "class  I*  occurrence. 


All  the  tests  on  this  module  are  limited  to  a  maximum  "N"  of  30,  which  is. 

somewhat  larger  than  Conover  used  in  exact  binomial  calculations.  To  use  a  larger 

"N",  you  may  either  redimension  the  arrays  {see  Instructions  on  Modifying  Program 

Modules  in  APPENDIX  C)  or  you  may  use  the  Normal  Approximation  recommended 

in  the  Conover  text  for  large  "N".  A  Normal  Approximation  is  not  included  in  the 

software  module.  The  tests  on  the  module  are  also  designed  to  accept  only  decimal 

values  for  "P",  the  hypothesized  probability,  (Note-  the  screen  format  on  the  MI00 
* 

uses  "P"  in  lieu  of  "P  ")  and  "ALPHA",  the  significance  level. 


If  you  wish  to  exit  a  program  in  progress,  strike  "shift"  and  keep  this  key  depressed 
while  striking  the  "break/pause"  button  above  the  keyboard,  or  the  "control"  key  and 
"C".  After  exiting  the  program  if  you  strike  "f4",  which  is  function  key  4,  or  type  "run" 
and  strike  the  "enter"  key  the  program  will  reinitialize.  If  you  wish  to  exit  to  the  Ml  00 
main  menu  after  breaking  out  of  the  program,  simply  strike  "f8",  which  is  function  key 
8,  or  type  the  word  "menu"  and  strike  the  "enter"  key. 


Once  you  enter  the  "BINOM"  module  the  Module  Main  Menu  will  appear  on  the 
screen  as  follows: 


1-Blnomlal  Test  2-Quantile  Test 

3-Cox-Stuart  Test  4-Quantile  C.  I. 

5-Siqn  Test 

PLEASE  MAKE  YOUk  SELECTION  <l-5> 


The  discussion  and  example  problems  for  the  tests  are  presented  in  this  user's  manual 
in  their  order  of  appearance  on  the  Binomial  Module  Main  Menu.  Each  test  is 
demonstrated  as  a  "self  contained  -  stand  alone"  example,  so  that  you  may  study  and 
use  them  in  any  order  you  like.  After  you  have  worked  through  several  of  the  tests 
within  the  module,  you  may  find  that  some  of  the  procedural  information  and  screen 
explanations  are  redundant. 
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!  Binomial  Test 


General 

The  Binomial  Test  is  used  to  determine  the  reasonableness  of  the  assumption 

* 

that  a  sample  could  have  been  drawn  from  a  population  having  the  value  "P  ",  which 
you  hypothesized  in  your  null  hypothesis. 

*  *  Example  Problem  *  * 

Problem  Setting 

Suppose  you  have  watched  people  rolling  a  certain  die  for  a  period  of  time,  and 
you  believe  that  the  "five"  is  coming  up  too  often.  You  wish  to  test  your  hypothesis 
statistically  to  determine  if  the  die  is  fair.  You  know  that  if  it  is,  then  all  the  faces  are 
equally  likely.  Your  null  hypothesis  for  any  face  should  be 

P*  =  P  =  1/6. 

"P  "  is  the  constant  probability  you  expect  if  the  null  hypothesis  is  true,  is  the 
probability  obtained  when  the  test  is  conducted,  and  1/6  is  the  probability  of  any  given 
face.  You  roll  the  die  "N"  trials  to  obtain  a  sample  to  use  to  test  your  hypothesis. 
Suppose  for  our  example  you  roll  the  die  12  times,  and  observe  that  the  "five"  occurred 
exactly  8  times.  Given  this  information  let  s  test  our  null  hypothesis. 

Computer  Solution  to  Problem 

One  possible  way  to  test  this  hypothesis  is  with  the  application  of  the  Binomial 
Test  to  the  data.  Therefore,  selection  "1"  is  made  from  the  Module  Main  Menu. 
(Note  -  it  is  not  necessary  to  strike  the  "enter"  key  after  making  your  selection.) 


1-Binomial  Test  2-Ouantile  Test 

3-Cox-Stuart  Test  4-Quantile  C.  I. 

5-Sign  Test 

PLEASE  MAKE  YOUR  SELECTION  <l-5>  1 


^  .s 
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The  following  screen  will  now  appear: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  1 Y 1  for  yes  or  1 N 1  for  no 

If  you  wish  to  check  the  assumptions,  simply  enter  a  Capital  "Y”  otherwise  you  must 
enter  a  Capital  "N".  ( The  machine  will  remain  in  the  loop  until  one  of  these  two  letters  is 
input.)  If  you  input  an  "N"  for  your  response,  the  M100  will  skip  the  data  assumption 
routine  and  request  your  input  parameters.  For  the  purposes  of  this  example  we  will 
assume  that  your  answer  was  "Y"  to  the  request  for  information  on  the  test 
assumptions.  The  following  screen  will  appear: 


You  may  strike  any  key  except  function  keys,  shift  keys,  or  control  keys  to 
continue. 

The  following  screen  now  will  appear: 


If  after  reading  the  assumptions,  you  decide  that  the  Binomial  would  not  be 
appropriate,  input  an  "M".  This  would  return  you  to  the  .Ylodule  Main  Menu  from 
which  you  may  select  another  test.  In  this  example  however,  we  find  that  the  Binomial 
lest  is  appropriate  to  perform  on  the  die  example  data  and  "strike  any  key"  to 


continue. 


The  program  is  now  ready  for  the  test  parameter  input.  The  first  parameter 
requested  is  "N". 


We  input  "12"  as  the  number  of  observations  or  trials  and  strike  the  "enter"  key.  After 

the  machine  has  accepted  the  first  parameter,  it  will  return  with  the  following  screen 

.  .  * 
display  requesting  that  you  input  "P  "  (the  hypothesized  P  value). 


We  input  the  decimal  value  for  our  hypothesized  probability  value  of  1/6. 

The  machine  now  requires  input  from  you  concerning  the  significance  level  or 
Alpha  level  at  which  you  wish  to  test  your  null  hypothesis.  The  program  will  accept  any 
decimal  values  between  0  and  I  for  Alpha.  Due  to  the  discrete  nature  of  the  Binomial 
Distribution  and  the  machine  program,  your  actual  alpha  will  always  be  "equal  to  or 
less  than"  the  alpha  value  which  you  input.  We  will  input  an  alpha  of  .075  for  our 
example. 


The  last  input  request  is  one  for  information  concerning  the  location  of  the  critical 
regions.  This  screen  requests  that  you  input  an  integer  indicating  where  you  want  the 
critical  region  located. 


3 


Please  select  your  critical  region 

1  -  Two  Tailed  Test 

2  -  Lower  Tal 1  Only 

3  -  Upper  Tall  Only 

PLEASE  MAKE  YOUR  SELECTION  <l-3>?  3 


We  elected  to  do  an  "Upper  Tailed  Test  ,  because  it  appeared  that  the  die  was  loaded  in 
favor  of  rolling  a  "five"  when  we  watched  it  being  thrown  before  we  conducted  our 
trials. 

Pressing  the  "enter"  key  displays  the  following  screen  which  is  self  explanatory. 


I'm  Calculating 


Upon  completion  of  the  calculations  the  screen  will  return  with  following  output 


ALPHA  is  equal  to  .036378 
reject  your  null  hypothesis  if 
' T '  is  greater  than  4 
TO  CONTINUE  -  STRIKE  ANY  KEY 


The  information  on  this  screen  is  the  output  data  from  the  Hypothesis  Test.  The  test 
statistic  T  is  defined  in  Conover  as  "the  number  of  times  the  outcome  is  ‘class  1" 
[Ref.  2:p.  97).  The  'T'  presented  here  indicates  the  number  of  “class  1"  outcomes  which 
lie  in  the  .036378  area  of  the  binomial  PDF.  After  striking  any  key,  the  M  100  will 
inquire  whether  you  wish  to  determine  what  your  ALPHA  HAT  value  is. 


Do  you  wish  to  find  your  actual 
ALPHA  HAT  VALUE? 

Just  Input  1 Y *  for  Yes  or  ‘ N 1  for  No 


The  MlOO  will  only  accept  a  Capital  T  or  a  Capital  "N".  Any  other  entry  will 
cause  the  computer  to  remain  in  the  loop  waiting.  If  the  computer  does  not 
immediately  respond,  check  to  see  that  you  have'Caps  Lock"  on.  DO  NOT  STRIKE 
THE  ENTRY  KEY  AFTER  INPUTING  THE  LETTER  OR  THE  COMPUTER  WILL 
SKIP  A  SCREEN  DISPLA  Y.  In  the  case  of  our  example  problem  we  will  input  a  "Y". 


Please  input  your  observed  1 T 1 ?  8 


You  are  asked  to  input  the  number  of  observed  "class  1"  occurrences.  In  our 
example  we  observed  8  trials  inwhich  "five"  appeared,  so  we  input  an  "8"  and  then 
press  the  "enter"  key. 

The  computer  responds  with  the  following: 


ALPHA  HAT  is  equal  to  1. 3E-05 
TO  CONTINUE  -STRIKE  ANY  KEY 


The  probability  of  having  "8"  or  more  "fives"  in  12  rolls  on  a  fair  die  is  a  very  rare 
event.  This  would  occur  by  chance  with  a  probability  of  only  .000013  if  the  null 
hypothesis  (i.e.  that  P=  1/6)  were  true.  When  we  reject  this  null  hypothesis  we  can  do 
so  with  relatively  strong  confidence  that  we  have  not  made  a  type  1  error. 
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We  have  now  completed  the  Binomial  Test,  and  strike  a  key.  The  following  screen 
appears: 


If  you  strike  "C",  the  machine  will  reset  itself.  Do  not  strike  the  "enter"  key  after 
making  your  selection. 

|  b.  1  Quantile  Test 

General 

The  quantile  test  is  an  application  of  the  binomial  test  to  an  ordinal  level,  discrete 
random  variable.  Examples  are  student's  test  scores  on  a  given  true/false  or  multiple 
choice  test.  You  could  use  the  Quantile  test  to  perform  an  hypothesis  test  on  sample 
data  taken  from  the  population  of  students  who  took  the  examination.  Using  the 
sample  data,  you  can  test  underlying  assumptions  which  you  may  have  concerning  the 
score  associated  with  a  certain  quantile  in  the  population  as  a  whole.  If  you  already 
know  what  scores  are  associated  with  the  quantiles  in  the  population,  you  may  wish 
instead  to  determine  if  the  sample  you  drew  from  the  population  is  representative  of 
the  whole  group.  If  the  quantile  associated  with  a  given  score  in  the  sample  is 
approximately  what  you  would  expect,  based  either  on  your  assumptions  or  upon  what 
you  know  about  the  group  scores  as  a  whole,  you  would  retain  your  hypothesis. 

*  *  Example  Problem  *  * 

Problem  Setting 

Suppose  you  have  data  on  Army  Physical  Fitness  Tests  (APFT)  from  three 
military  posts.  You  know  the  combined  median  scores  on  each  of  the  three  events 
were  15.6  minutes  for  the  2  mile  run,  43  pushups,  and  51  situps. 


You  have  just  tested  a  unit  of  10  soldiers  and  their  scores  on  the  pushups  listed  below 
appear  to  be  higher  than  the  "Army  Average": 

38  43  43  46  51  59  60  60  61  62 

Since  this  is  a  Special  Forces  unit,  you  would  expect  the  median  score  achieved  in  this 
unit  to  be  higher  than  the  "Army  Average".  Formally  stated  our  hypothesis  is  as 
follows: 


Hq.  "SF"  median  ^  "AA"  median 
Hj.  "SF"  median  >  "AA"  median 

Computer  Solution  to  the  Problem 

An  examination  of  the  data  would  tend  to  indicate  that  the  Quantile  Test  might  be 
appropriate.  So  we  chose  selection  "2"  from  the  Module  Main  Menu. 


MENU 

1-Binomial  Test  2-Ouantile  Test 

3-Cox-Stuart  Test  4-Quantile  C. I. 

5-Sign  Test 

PLEASE  MAKE  YOUR  SELECTION  <l-5>  2 


(Note  that  it  is  not  necessary  to  strike  the  "enter"  key  after  selecting  the  test  you  wish  to 
perform  .)  The  following  screen  appears: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  * Y 1  for  yes  or  1 N 1  for  no 


For  the  example,  we  will  check  the  assumptions  of  the  Quantile  l  est  by  inputing  a 
Capital  "Y".  Inputing  a  Capital  "N"  causes  the  M100  to  skip  the  data  assumption 
routine  and  request  your  input  parameters.  This  will  provide  you  with  the  following 
two  screens: 
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If  you  fail  any  of  the  assumptions,  then  you  should  input  an  M  to  return  to  the 
Module  Main  Menu.  In  this  example,  we  meet  the  requirements  for  a  Quantile  Test. 
So  we  strike  a  key  proceed  with  our  test. 

The  program  now  requests  that  you  input  parameters  to  be  used  in  the  hypothesis 
test  to  be  performed.  The  first  parameter  requested  for  input  is  N. 


Please  input  N?  10 


Ten  people  took  the  PT  test,  so  we  input  the  number  "10"  and  strike  the  "enter"  key. 

♦ 

The  next  screen  requests  that  you  input  P  ,the  hypothesized  probability.  In  the 
example  problem,  we  are  interested  in  testing  the  "Median",  so  we  would  input  ".5". 


M  t.1  Ja»  Jy A-  !»«.>,  b*. 


The  M100  will  now  ask  you  for  information  concerning  the  significance  level  or 
Alpha  Level  The  program  will  accept  any  decimal  values  between  0  and  1  for  Alpha. 
Due  to  the  discrete  nature  of  the  Binomial  Distribution  and  the  machine  program,  your 
actual  alpha  will  be  "equal  to  or  less  than"  the  alpha  value  which  you  input.  We  will 
arbitrarily  select  an  Alpha  of  .10  for  this  example. 


Please  Input  the  ALPHA  VA1UE  .10 


The  final  request  for  parameters  concerns  the  location  of  the  critical  regions.  You 
can  chose  a  two  tailed  or  a  one  tailed  test.  This  step  should  be  performed  carefully, 
since  it  is  easy  to  make  an  error  here  when  doing  a  Quantile  Test.  We  believe  that  our 
small  group  of  ten  individuals  will  have  a  median  score  higher  than  the  "Army 
Average",  so  we  need  to  perform  a  one  tailed  hypothesis  test.  We  expect  that  if  we 
rank  order  the  scores  of  our  sample  and  calculate  a  CDF,  the  percentile  for  the  Army 
Average  median  score  (43)  will  be  lower  than  .5  when  that  score  is  encountered  in  our 
sample  set.  Therefore  we  select  a  "Lower  Tailed  Test"  and  enter  selection  "2"  as  our 
choice. 


Please  select  your  critical  region 

1  -  Two  Tailed  Test 

2  -  Lower  Tail  Only 

3  -  Upper  Tai 1  Only 


PLEASE  MAKE  YOUR  SELECTION  <l-3>?  2 


We  have  made  our  selection  now  and  strike  the  "enter"  kev. 


While  the  Binomial  Probability  computation  is  in  progress  the  screen  will  display: 


N  \  >  * 


i 

J 


Upon  completion  of  the  calculations  the  following  screen  will  appear: 


ALPHA  Is  equal  to  .0546875 
reject  your, null  hypothesis 


'TVis  le'ss  than  or  e 
TO  CONTINUE  -  STRIKE 


ual  to 
KEY 


if 

2 


The  "'Tl'"  is  the  number  of  observations  which  are  "less  than  or  equal  to  X  [Ref.  2:p. 
107].  This  screen  would  be  different  if  you  had  selected  an  upper  tailed  or  a  two  tailed 
test.  In  the  case  of  an  upper  tailed  test,  "'T2'"  would  be  displayed  rather  than  '"Tl"'. 

rjc" 

The  "T2"'  value  is  the  number  of  observations  that  are  "less  than  X  [Ref.  2:p.  107], 
You  would  reject  your  null  hypothesis  if  your  observed  "'T2'"  value  was  larger  than  the 
screen  value  "'T2'"  displayed.  If  a  two  tailed  test  had  been  selected,  both  "'Tl'"  and 
"'T2'"  would  have  been  displayed.  You  would  be  seeking  either  values  smaller  than 
"TT"  or  larger  than  After  striking  any  key,  the  M100  will  inquire  whether  you 

wish  to  determine  what  you  ALPHA  HAT  value  is. 


Do  you  wish  to  find  your  actual 
ALPHA  HAT  VALUE? 

Just  input  ‘Y1  for  Yes  or  'N1  for  No 


If  your  response  were  Capital  ",Y"  then  the  M1(H)  would  simply  go  to  the  final 
screen.  For  the  example,  we  will  request  the  actual  APLHA  HAT  value,  so  we  input  a 
Capital  Y.  (Note  -  you  do  not  need  to  strike  the  "enter"  key  as  the  machine  will 

accept  the  input  without  it.)  The  machine  requires  the  input  of  your  actual  observed 

* 

"Tl  "  and  ”T2'".  In  this  example  X  is  the  "Army  Average"  median  score.  In  our 
data  on  the  number  of  pushups,  we  had  "3"  scores  equal  to  or  less  than  the  median 
score  of  43,  and  "\"  score  less  han  43.  These  values  are  input  on  two  successive 
screens  which  appear  as  follows: 


Please  input  your  observed  'Tl?  3 


Please  input  your  observed  ' T2 1 ?  1 


(Note  -  inputing  a  "'7 2'"  value  which  is  larger  than  "'77'"  gives  a  machine  error.) 

The  computer  responds  with  the  following  output  screen  indicating  ALPHA  HAT, 
based  on  the  observed  data.  The  probability  of  having  a  "'Tl'"  value  less  than  or  equal 
to  3,  if  the  true  median  score  is  43,  would  occur  by  chance  with  probability  of  .171875. 
Our  "'Tl'"  value  was  too  large  to  reject  our  null  hypothesis  at  an  alpha  level  smaller 
than  .171875. 


We  have  completed  the  Quantile  Test.  We  strike  a  key  and  the  following  screen 
appears: 


Strike  'C1  to  do  another  test 
If  you  are  finished  strike  any  other  key 


If  you  strike  "C",  the  machine  resets  itself  to  the  Module  Main  Menu.  The  input  of  any 
other  key  prepares  the  machine  to  receive  either  the  typed  entry  "Menu",  or  the  input 
of  "f8"  function  key  to  return  to  the  M100  main  menu. 


|Sign  Test 


General 


The  Sign  Test,  as  treated  by  Conover,  is  based  on  a  binomial  test  with  P  set  to 
1/2.  The  sign  test  is  used  with  related  samples.  The  most  common  of  these  related 
samples  are  either  "before"  and  "after"  data  samples  or  matched  pairs.  Let's  examine  a 
problem  which  will  explain  the  sign  test  procedures. 

*  *  Example  Problem  *  * 

Problem  Setting 

Say  you  have  10  sets  of  identical  twins,  and  you  believe  that  "first  born"  twins  are 
more  aggressive  and  dominant  than  the  "second  born".  To  test  this  you  develop  and 
administer  a  psychological  profile  battery  to  get  an  aggression  scale  for  each  of  the  20 
people  in  the  survey  After  completing  your  analysis,  you  got  the  following  scores  on  a 
scale  from  1  to  10  on  the  aggression  index. 


(First  Born) 
Y  SCORE 


(Second  Born)  Associated 
X  SCORE  Si qn 


You  are  interested  in  determining  if  our  inital  hypothesis  is  true.  Stated  formally, 
the  hypothesis  would  be  as  follows: 

H0;  P(  +  )  <  P(-) 

H1:  P(  +  )  >  P(-) 


Computer  Solution  to  Problem 


We  have  ordinal  level  data  from  10  matched  pairs;  therefore,  we  select  "5"  from 
the  Module  Main  Menu.  (Note  that  it  is  not  necessary  to  strike  the  "enter"  key  after 
selecting  your  test.) 


MENU 

1-Binomial  Test  2-Ouantile  Test 

3-Cox-Stuart  Test  4-Quantile  C. I. 

5-Sign  Test 

PLEASE  MAKE  YOUR  SELECTION  <l-5>  5 


The  M 100  now  asks  you  if  you  wish  to  check  our  assumptions. 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  'Y1  for  yes  or  1 N '  for  no 


If  you  wish  to  check  the  assumptions,  simply  enter  a  Capital  "Y"  otherwise  you  must 
enter  a  Capital  "AT.  ( The  machine  will  remain  in  the  loop  until  one  of  these  two  letters  is 
input.)  If  you  input  an  "N"  for  your  response,  the  M100  will  skip  the  data  assumption 
routine  and  request  your  input  parameters.  For  the  purposes  of  this  example  we  will 
assume  that  our  answer  was  "Y"  to  the  request  for  information  on  the  test 
assumptions.  The  following  screen  will  appear: 


You  selected  the  Sign  Test 
this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


At  this  time  you  may  strike  any  key  except  the  function  keys  located  above  the 
keyboard,  the  shift  keys,  or  the  control  keys  in  order  to  continue. 


The  following  screen  listing  the  assumptions  of  the  sign  test  will  now  appear. 


We  input  a  "1"  here  because  there  were  seven  times  that  the  "X"  aggression 
scores  were  lower  than  the  "Y"  aggression  score. 

The  next  screen  appears  which  requests  the  information  concerning  the  reverse 
situation. 


We  input  an  "\"  and  strike  the  "enter"  key  to  continue.  We  input  this  number 
because  among  our  10  sets  of  twins,  in  only  one  case  did  the  second  born  twin  have  a 
higher  aggression  score. 

Now  the  M100  will  request  input  from  you  concerning  the  significance  level  or 
Alpha  Level  at  which  you  wish  to  test  your  null  hypothesis.  It  will  accept  any  decimal 
value  between  0  and  1  for  alpha. 


Please  input  the  ALPHA  VA1UE  .10 


The  Binomial  Distribution  is  discrete,  and  the  machine  has  been  programmed  to  ensure 
that  your  actual  alpha  will  always  be  "equal  to  or  less  than"  the  alpha  value  which  you 
input.  We  input  .10  for  our  example. 

The  final  input  request  concerns  the  location  of  your  critical  regions.  This  screen 
requests  you  input  an  integer  indicating  where  you  want  the  critical  region  located. 
You  can  chose  a  two  tailed  test  or  a  one  tailed  test.  This  step  should  be  performed 
carefully,  as  the  assignment  of  X  and  Y  directly  influence  the  location  of  the  critical 
region.  The  tail  that  we  wish  to  examine  is  the  "Upper  Tailed  Test"  because  we  want 
to  place  our  original  conjecture  on  trial. 


Please  select  your  critical  region 

1  -  Two  Tailed  Test 

2  -  Lower  Tail  Only 

3  -  Upper  Tail  Only 

PLEASE  MAKE  YOUR  SELECTION  <l-3>?  3 


The  next  screen  is  self  explanatory. 


I'm  Calculating 


Upon  completion  of  the  calculations  the  following  screen  will  appear: 


ALPHA  is  equal  to  .0351563 
reject  your  null  hypothesis 
T  is  greater  than  or  equal  to 
TO  CONTINUE  -  STRIKE  ANY  KEY 


if 

7 


This  screen  would  have  appeared  slightly  different  if  we  had  chosen  a  lower  tailed  test 
or  a  two  tailed  test.  After  striking  any  key,  The  M 100  will  inquire  whether  you  wish  to 
determine  what  your  ALPHA  HAT  value  is. 


Do  you  wish  to  find  your  actual 
ALPHA  HAT  VALUE? 

Just  input  1 Y 1  for  Yes  or  1 N 1  for  No 


The  M100  will  only  accept  a  Capital  "  Y"  or  a  Capital  "iX".  Any  other  entry  will 
cause  the  computer  to  remain  in  the  loop  waiting.  If  the  computer  does  not 
immediately  respond,  check  to  see  that  you  have  "Caps  Lock"  on.  DO  SOT  STRIKE 
THE  ESTRY  KEY  AFTER  LXPETIXG  THE  LETTER  OR  THE  COMPUTER  WILL 
SKIP  A  SCREES  DISPLA  Y.  In  the  case  of  our  example  problem  we  will  input  a  "Y". 


Please  input  your  observed  1 T 1 ?  7 


The  value  to  input  for  "T"'  is  the  number  of  observed  "Pluses".  We  observed  "7" 
cases,  so  this  is  the  number  we  input. 

The  computer  responds  with  the  following. 


ALPHA  HAT  is  equal  to  .0351563 
TO  CONTINUE  -STRIKE  ANY  KEY 


This  means  that  by  chance,  seven  or  more  Pluses  would  occur  with  a  probability  of 
.0351563  under  a  true  null  hypothesis. 

We  have  completed  the  Sign  Test/Cox-Stuart  Test.  Strike  any  key  and  the 
following  screen  appears. 


Strike  'C1  to  do  another  test 
If  you  are  finished  strike  any  other  key 


If  you  strike  "C",  the  machine  resets  itself  to  the  Module  Main  Menu.  The  input  of  any 
other  key  prepares  the  machine  to  receive  either  the  typed  entry  "Menu",  or  the  input 
of  "f8"  function  key  to  return  to  the  M100  main  menu. 

|d.  ICox-Stuart  Test 

General 

The  Cox-Stuart  test  is  very  similar  to  the  sign  test.  It  also  is  based  on  a  binomial 
test  with  a  P  equal  to  1/2.  The  major  difference  lies  in  the  fact  that  the  Cox-Stuart  test 
is  used  either  as  a  test  for  trends,  or  as  an  ordinal  level  test  of  correlation  on  a  single 
sample. 


*  *  Example  Problem  *  * 


Problem  Setting 

Let  us  assume  that  you  have  a  sample  of  20  people  of  various  ages.  We  believe 
that  as  people  become  older,  they  tend  to  become  less  flexible,  or  more  set  in  their 
ways.  To  test  this  you  develop  and  administer  a  psychological  profile  battery  to 
determine  the  degree  of  flexiblity  of  these  20  people.  The  scores  you  derived  are  as 
follows: 


CASE 

AGE 

SCORE 

1 

26 

3 

2 

27 

6 

3 

34 

3 

4 

35 

2 

5 

39 

4 

6 

41 

2 

7 

42 

7 

8 

42 

2 

9 

45 

9 

10 

53 

4 

11 

54 

6 

12 

54 

5 

13 

57 

8 

14 

58 

7 

15 

61 

6 

16 

61 

5 

17 

65 

3 

18 

71 

9 

19 

73 

6 

20 

74 

2 

We  divide  these  sample  scores  into  two  halves.  In  this  manner,  we  form 
(Xj,Xj  +  c)  pairs.  These  are  now  treated  the  same  as  the  (X,Y)  pairs  were  in  the  sign 
test.  A  "Plus"  is  assigned  when  the  X  score  is  lower  than  the  Y  score,  and  a  "Minus"  is 
assigned  when  the  X  score  is  higher  than  the  Y  score  in  the  pair. 

(Younger)  (Older)  Associated 

PAIR  X  SCORE  Y  SCORE  Sian 


You  are  interested  in  determining  if  our  inital  hypothesis  is  true.  Stated  formally, 
the  hypothesis  would  be  as  follows: 

H0.  P(  +  )  ^P(-) 

H,;  P(  +  )  >  P(-) 

Computer  Solution  to  Problem 

We  have  ordinal  level  data  and  we  are  looking  for  a  trend;  therefore  we  select  "3" 
from  the  Module  Main  Menu.  (Note  that  it  is  not  necessary  to  strike  the  "enter"  key 
after  selecting  your  test.) 

MENU 

1-Binomial  Test  2-Quantile  Test 

3-Cox-Stuart  Test  4-Quantile  C. I. 

5-Sign  Test 

PLEASE  MAKE  YOUR  SELECTION  <l-5>  3 

The  M 100  now  asks  you  if  you  wish  to  check  our  assumptions. 

Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  ' Y 1  for  yes  or  1 N 1  for  no 

If  you  wish  to  check  the  assumptions,  simply  enter  a  Capital  "Y"  otherwise  you  must 
enter  a  Capital  " N".  ( The  machine  will  remain  in  the  loop  until  one  of  these  two  letters  is 
input.)  If  you  input  an  "N"  for  your  response,  the  MI00  will  skip  the  data  assumption 
routine  and  request  your  input  parameters.  For  the  purposes  of  this  example  we  will 
assume  that  our  answer  was  "Y"  to  the  request  for  information  on  the  test 
assumptions.  The  following  screen  will  appear: 

You  selected  the  Cox-Stuart  Test 
this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


At  this  time  you  may  strike  any  key  except  the  function  keys  located  above  the 
keyboard,  the  shift  keys,  or  the  control  keys  in  order  to  continue. 

The  following  screen  listing  the  assumptions  of  the  sign  test  will  now  appear. 

1- Ordinal  or  Higher  Level  Data 

2- All  Observations  are  I ID 

TO  CONTINUE  -  STRIKE  ANY  KEY 
If  you  fail  the  assumptions  then 
to  return  to  the  menu  input 

The  flexibility  scale  you  developed  is  certainly  ordinal,  and  the  data  points  are 
independent  and  identically  distributed.  It  would  therefore  appear  that  we  meet  the 
assumptions  of  the  test.  .You  continue  by  "striking  a  key". 

The  machine  will  now  query  you  for  the  parameters  needed  to  perform  the  test.  In 
the  sign  test,  you  are  interested  in  the  number  of  times  "X  is  Less  than  Y"  and  "X  is 
Greater  than  Y" .  In  this  example,  we  classified  the  younger  person  as  X,  and  the  older 
as  Y.  This  classification  was  arbitrary.  If  the  classification  were  reversed,  the  choice  of 
a  tail  location  on  a  one  tailed  test  would  have  to  be  switched.  If  we  had  had  ties  where 
"X  is  equal  to  Y"  these  would  have  been  omitted. 


Input  the  #  of  times  X<Y  (+)  6 


We  input  a  "6"  here  because  there  were  six  cases  that  younger  person's  flexibility 
scores  were  lower  than  the  older  person's  score. 


UPJWJ.U 
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The  next  screen  appears  which  requests  the  information  concerning  the  reverse 
situation. 


We  input  an  "4"  and  strike  the  "enter"  key  to  continue. 

Now  the  M100  will  request  input  from  you  concerning  the  significance  level  or 
Alpha  Level  at  which  you  wish  to  test  your  null  hypothesis.  It  will  accept  any  decimal 
value  between  0  and  1  for- alpha. 


The  Binomial  Distribution  is  discrete,  and  the  machine  has  been  programmed  to  ensure 
that  your  actual  alpha  will  always  be  "equal  to  or  less  than"  the  alpha  value  which  you 
input.  We  input  .10  for  our  example. 

The  final  input  request  concerns  the  location  of  your  critical  regions.  This  screen 
requests  you  input  an  integer  indicating  where  you  want  the  critical  region  located. 
You  can  chose  a  two  tailed  test  or  a  one  tailed  test.  As  with  the  Sign  test  discussed 
above,  this  step  should  be  performed  carefully,  as  the  assignment  of  X  and  Y  directly 
influence  the  location  of  the  critical  region.  The  tail  that  we  wish  to  examine  is  the 
"Upper  Tailed  Test"  because  we  want  to  place  our  original  conjecture  on  trial. 


Please  select  your  critical  region 

1  -  Two  Tailed  Test 

2  -  Lower  Tai 1  Only 

3  -  Upper  Tail  Only 

PLEASE  MAKE  YOUR  SELECTION  <l-3>?  3 
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We  input  a  "3"  for  the  upper  tailed  test. 
The  next  screen  is  self  explanatory. 


I'm  Calculatl ng 


Upon  completion  of  the  calculations  the  following  screen  will  appear: 


ALPHA  is  equal  to  .0546875 
reject  your  null  hypothesis  i 
1 T  is  greater  than  or  equal  to 
TO  CONTINUE  -  STRIKE  ANY  KEY 


f 
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This  screen  would  have  appeared  slightly  different  if  we  had  chosen  a  lower  tailed  test 
or  a  two  tailed  test.  After  striking  any  key,  The  M100  will  inquire  whether  you  wish  to 
determine  what  your  ALPHA  HAT  value  is. 


Do  you  wish  to  find  your  actual 
ALPHA  HAT  VALUE? 

Just  input  1 Y '  for  Yes  or  1 N 1  for  No 


The  MlOO  will  only  accept  a  Capital  "Y"  or  a  Capital  "N".  Any  other  entry  will 
cause  the  computer  to  remain  in  the  loop  waiting.  If  the  computer  doc  not 
immediately  respond,  check  to  see  that  you  have  "Caps  Lock"  on.  DO  NOT  STRIKE 
THE  ENTRY  KEY  AETER  INPUTING  THE  LETTER  OR  THE  COMPUTER  WILL 
SKIP  A  SCREEN  DISPLA  Y.  In  the  case  of  our  example  problem  we  will  input  a  "Y". 


The  value  to  input  for  '"T”'  is  the  number  of  observed  "Pluses".  We  observed  "6" 
cases,  so  this  is  the  number  we  input. 

The  computer  responds  with  the  following. 


ALPHA  HAT  is  equal  to  .3769533 
TO  CONTINUE  -  STRIKE  ANY  KEY 


This  means  that  by  chance,  six  or  more  Pluses  would  occur  with  a  probability  of 
.376953  under  a  true  null  hypothesis.  This  outcome  does  not  indicate  a  trend  to 
become  less  flexible  as  one  becomes  older,  so  we  would  retain  our  null  hypothesis. 

We  have  completed  the  Cox-Stuart  Test.  Strike  anv  key  and  the  follov  ;ng  screen 
appears. 


Strike  1 C 1  to  do  another  test 
If  you  are  finished  strike  any  other  key 


If  you  strike  "C",  the  machine  resets  itself  to  the  Module  Main  Menu.  The  input  of  any 
other  key  prepares  the  machine  to  receive  either  the  typed  entry  "Menu",  or  the  input 
of  “f8"  function  key  to  return  to  the  M100  main  menu. 


I  e.  |  Quantile  Confidence  Interval 

General 

The  Quantile  Confidence  Interval  test  is  an  application  of  the  binomial  test  to  an 
ordinal  level,  discrete  random  variable.  The  test  permits  you  to  determine  a  confidence 
interval  equal  "7  -  Alpha". 

*  *  Example  Problem  *  * 

Problem  Setting 

Recall  that,  in  the  Quantile  Test  example,  we  had  data  on  Army  Physical  Fitness 
Tests  from  three  military  posts.  We  were  given  the  combined  median  scores  on  each  of 
the  three  events.  These  median  scores  were  15.6  minutes  for  the  2  mile  run,  43 
pushups,  and  51  situps.  We  tested  a  Special  Forces  unit's  pushup  scores  to  determine 
if  they  were  really  better  than  the  "Army  Average"  soldier.  Now  you  have  been 
directed  by  your  commander  to  examine  the  scores  of  one  specific  28  man  platoon. 
This  platoon  has  been  identified  as  seeming  to  have  substandard  performance.  The 
commander  tells  you  that  the  statistical  findings  will  be  sent  to  higher  headquarters  and 
will  be  the  determining  factor  in  whether  or  not  the  noncommissioned  officer  in  charge 
is  permitted  to  remain  in  the  Army.  The  commander  would  therefore  like  to  be  at  least 
95%  confident  that  this  is  not  a  chance  occurrence.  The  scores  on  the  pushups  for  this 
platoon  are  as  follows: 

34  36  36  37  37  37  38  38  39  39  39  40  41  41 
42  43  43  43  45  45  46  48  49  52  58  60  60  62 
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Computer  Solution  to  the  Problem 

It  seems  appropriate,  based  on  the  directions  from  your  boss,  to  calculate  a 
Quantile  Confidence  Interval.  So  we  chose  selection  "4"  from  the  Module  Main  Menu. 


1-Binomial  Test  2-0uantile  Test 

3-Cox-Stuart  Test  4-Quantile  C.  I. 

5-Sign  Test 

PLEASE  MAKE  YOUR  SELECTION  <l-5>  4 


( Note  that  it  is  not  necessary  to  strike  the  "enter"  key  after  selecting  the  test  you  wish  to 
perform.)  The  following  screen  appears: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  1 Y '  for  yes  or  1 N 1  for  no 

For  the  example,  we  will  check  the  assumptions  of  the  Quantile  Confidence  Interval, 
by  inputing  a  Capital  "  Y".  If  you  have  already  performed  the  Quantile  Test  you  would 
probably  input  a  Capital  "N",  because  the  assumptions  are  identical  for  both  tests,  and 
you  are  using  the  same  type  of  data.  This  "N"  input  would  causes  the  M100  to  skip 
the  data  assumption  routine  and  request  your  input  parameters.  If  you  input  "Y"  the 
M100  provides  the  following  two  screens: 


You  selected  the  Quantile  Test/CI 
this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


1- Ordinal  or  Higher  Level  Data 

2- All  Observations  Are  I  ID 

TO  CONTINUE  -  STRIKE  ANY  KEY 
If  you  fail  the  assumptions  then 
to  return  to  the  menu  input  Mr 


If  you  fail  any  of  the  assumptions,  then  you  should  input  an  M  to  return  to  the  Main 
Menu.  In  this  example,  we  meet  the  requirements  for  a  Quantile  Confidence  Interval 
Test.  So  we  strike  a  key  proceed  with  our  test. 

The  program  now  requests  the  four  parameters  necessary  to  perform  the 
Confidence  Interval  calculations.  The  first  is  N. 


Please  input  N?  28 


We  know  we  are  interested  in  a  28  man  platoon  so  we  input  that  number. 

The  next  screen  requests  that  you  input  a  decimal  value  between  0  and  1  for  the 
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P  (hypothesized  probability).  We  were  only  given  information  on  the  "Median" 
scores,  so  we  will  use  ".5  "  for  our  probability. 


Please  input  P(robabi 1 ity)  .50 


We  strike  the  "enter"  key  and  continue. 

The  M100  will  now  ask  you  for  information  concerning  the  significance  level  or 
Alpha  Level.  The  program  will  accept  any  decimal  values  between  0  and  1  for  Alpha. 
Due  to  the  discrete  nature  of  the  Binomial  Distribution  and  the  machine  program,  your 
actual  alpha  will  be  "equal  to  or  less  than"  the  alpha  value  which  you  input.  We 
arbitrarily  select  an  alpha  of  .05  for  this  example.  This  gives  us  a  ^  95%  confidence 
interval. 


The  next  screen  presents  you  with  the  actual  CDF  values  for  three  observations. 
The  observation  which  is  closest  to  the  requested  alpha  value  you  input  earlier  and  and 
its'  respective  CDF  value  is  presented  as  the  middle  case.  Here  we  had  input  alpha 
equal  to  ".05",  and  the  closest  alpha  value  in  the  discrete  binomial  distribution  occurs 
at  the  9th  observation  and  has  a  value  of  .0435793.  The  next  smaller  and  the  next 
larger  CDF  values  are  also  presented.  We  chose  this  method  of  presentation  to  permit 
you  to  decide  if  you  wished  to  use  an  alpha  value  which  was  slightly  higher  or  slightly 
larger  than  that  which  you  had  originally  selected.  In  this  case,  .0435793  is  close  to 
.05,  so  we  select  "9"  as  our  N  value,  and  strike  the  "enter"  key. 


LOWER  CDF 
. 0178491 
.0435793 
.0924667 


Input  the  N  for  the  lower  0BS#  9 


If  you  had  selected  an  upper  tailed  confidence  interval,  the  observation  numbers  would 
have  been  associated  with  their  respective  "  1  -  CDF'  value  and  you  would  have  selected 
an  "upper"  observation  number.  The  two  tailed  confidence  interval  presents  both  the 
upper  and  lower  observations  associated  with  their  "CDF'  and  "I  -  CDF'  values.  In 
that  case,  you  would  select  an  upper  and  a  lower  observation  number.  Your  actual 
alpha  value  would  be  the  sum  of  the  two  associated  values. 

The  next  screen  appears  with  the  data  we  need  to  assist  the  commander. 


ALPHA  is  equal  to  .0435793 

1 T1 1  lower  C. I.  is  <  or  =  10 
TO  CONTINUE  -  STRIKE  ANY  KEY 


T  his  is  the  information  you  need.  This  screen  indicates  that  if  you  have  10  or  fewer 

❖  * 

scores  that  are  "equal  to  or  less  than  X  ",  where  X  is  the  hypothesized  median,  then 
you  cannot  reject  that  X+  may  be  the  true  median.  In  our  example,  there  arc  18  scores 
which  are  less  than  or  equal  to  "43".  T  his  too  extreme  a  number  to  have  occurred  by 


chance,  if  the  median  of  the  sample  data  was  really  "43".  You  should  inform  the  boss 
that  this  platoon  is  does  not  appear  to  be  as  physically  fit  as  the  "Army  Average",  and 
that  he  can  be  more  than  .95°/o  confident  that  the  median  of  the  sample  is  less  than 
"43".  (If  this  had  been  a  two  tailed  interval,  you  would  have  had  a  screen  appear 
which  asked  if  you  wanted  the  machine  to  determine  the  value  of  the  data  points 
associated  with  the  "'Tl'"  and  the  "'T2'"  values.  An  answer  of  yes  means  that  the  M 100 
will  ask  for  you  to  input  all  your  scores.) 

We  have  completed  the  Quantile  Confidence  Interval  Test.  We  strike  a  key  and 
the  following  screen  appears: 


Strike  ‘C  to  do  another  test 
If  you  are  finished  strike  any  otner  key 


If  you  strike  "C",  the  machine  resets  itself  to  the  Module  Main  Menu.  The  input  of  any 
other  key  prepares  the  machine  to  receive  cither  the  typed  entry  "Menu",  or  the  input 
of  "IS"  function  key  to  return  to  the  M100  main  menu. 


3.  TESTS  BASED  ON  THE  CHI  SQUARE  DISTRIBUTION 

In  the  "CI1ISQ"  module  are  the  Cochran  test,  median  test,  chi  square  test, 
McNcmar  (  N  ^  20),  chi  square  goodness  of  fit  test,  and  tolerance  limits.  These  tests 
arc  all  derived  from  the  chi  square  distribution,  which  is  a  sum  of  "N"  independent 
standard  normal  variables,  squared  (Ref.  6].  The  degrees  of  freedom  associated  with 
the  chi  square  test,  are  a  function  of  the  number  of  squared  variables.  When  the 
degrees  of  freedom  equal  one,  the  square  root  of  the  chi  square  statistic  is  identical  to 
the  Z  score  for  a  normal(O.l). 

This  module  is  designed  so  that  when  ALPHA  is  requested,  it  must  be  input  as  a 
decimal  value.  The  CIIISQ  module  also  limits  the  Chi  Square,  McNemar,  and  the 
Median  test  to  contingency  tables  that  are  no  greater  than  14  rows  by  14  columns. 
To  use  a  larger  table,  see  APPENDIX  C  of  this  paper  for  instructions  on 
redimensioning  the  arrays. 
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If  you  wish  to  exit  this  program  when  in  progress,  strike  "shift"  and  keep  this  key 
depressed  while  striking  the  "break/pause"  button  above  the  keyboard,  or  strike  the 
"control"  key  and  "C".  If  you  wish  at  that  point  to  exit  to  the  M100  main  menu,  you 
need  only  strike  the  "f8"  function  key.  If  you  strike  the  "f4"  function  key,  the  program 
will  reinitialize  and  run  another  test. 

Once  you  enter  the  "CHISQ"  module  and  the  program  has  initialized,  the  Module 
Main  Menu  will  appear  on  the  screen  as  follows: 


MENU 

1-Tolerance  Limits 

2-McNemar  Test 

3-Cochran  Test 

4-Median  Test 

5-Chi  Square  Test 

6-Goodness  of  fit 

PLEASE  MAKE  YOUR  SELECTION  <l-6> 

Each  procedure  in  the  module  is  demonstrated  as  a  "self  contained  -  stand  alone" 
example,  so  that  you  may  study  and  use  them  in  any  order  you  like.  After  you  have 
worked  through  several  of  the  tests  within  the  module,  you  may  find  that  some  of  the 
procedural  information  and  screen  explanations  are  redundant. 

|  a.  |  Tolerance  Limits 

General 

Tolerance  limits  provide  an  interval  within  which  we  can  be  at  least  "1-alpha" 
confident  that  at  least  a  proportion  of  the  population  equal  to  "q"  is  within  the  interval 
[Ref.  2:p.ll7).  One  common  use  of  tolerance  limits  is  to  determine  the  sample  "N" 
needed  to  ensure  that  the  above  criteria  are  met.  Although  the  test  makes  no 
assumptions  about  the  underlying  distribution,  it  docs  assume  cither  that  there  is  an 
infinite  population  from  which  the  sample  is  drawn,  or  that  sampling  is  done  with 
replacement. 
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*  *  Example  Problem  *  * 


Problem  Setting 

Suppose  you  work  in  a  section  that  tests  expendable  military  equipment.  Your  job 
requires  that  you  determine  how  large  an  "N"  should  be  selected  out  of  equipment  lots 
in  order  to  be  confident  that  the  manufacturer  is  providing  equipment  which  meets  the 
contract  specifications.  The  most  recent  item  with  which  you  are  concerned  is  that  of 
nylon  helicopter  cargo  slings.  The  contractor  guarantees  that  with  95%  reliability 
these  slings  will  meet  their  specified  load  rating.  Your  boss  wants  to  be  "confident" 
that  this  is  so.  You  determine  that  this  means  your  boss,  who  is  not  a  trained 
operations  analyst,  would  probably  accept  an  "ALPHA"  value  or  Type  I  error  of  ".05". 

.Computer  Solution  to  the  Problem 

We  wish  to  establish  the  size  of  the  sample  required  to  ensure  that  at  least  a 
certain  percentage  lies  within  a  certain  confidence  interval.  Therefore  we  will  select  "l  " 
from  the  Module  Main  Menu.  ( Note  -  it  is  not  necessary  to  strike  the  "enter"  key  after 
making  your  selection.) 

MENU 

1-Tolerance  Limits  2-McNemar  Test 
3-Cochran  Test  4-Median  Test 

5-Chi  Square  Test  6-Goodness  of  fit 

PLEASE  MAKE  YOUR  SELECTION  <l-6>  1 


The  following  screen  will  now  appear: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  ’ Y 1  for  yes  or  1 N 1  for  no 


If  you  wish  to  check  the  assumptions,  simply  enter  a  Capital  "Y"  otherwise  a  Capital 
A  cause  the  machine  to  skip  the  data  assumption  loop.  Be  certain  that  you  have 
<  IPS  link  .  n  -alien  you  select  the  letter  otherwise  the  M100  will  remain  in  the 
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loop  indefinitely.  We  will  assume  that  you  answer  an  "Y"  to  the  request  for 
information  on  the  test  assumptions.  The  following  screen  will  appear: 


You  selected  the  Tolerance  Limits 
this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


You  may  strike  any  key  except  function  keys,  shift  keys,  or  control  keys  to 
continue. 


The  following  screen  now  will  appear: 


1- The  X's  are  a  random  sample 

2- Ordinal  level  data 

TO  CONTINUE  -  STRIKE  ANY  KEY 
If  you  fail  the  assumptions  then 
to  return  to  the  menu  input  M‘ 


If  you  decide  that  you  do  not  meet  the  assumptions  of  the  test  then  input  an  " M" .  This 
returns  you  to  the  Module  Main  Menu,  from  which  you  may  select  another  test.  For 
our  example,  we  decide  that  we  have  chosen  the  correct  test,  and  "strike  any  key"  to 
continue. 

The  program  is  now  ready  for  the  test  parameters.  To  calculate  tolerance  limits, 
we  must  input  four  parameters.  These  are  designated  as  "q",  "r" ,  "m" ,  and  "I  -  alpha". 
1  he  first  parameter  requested  is  the  confidence  interval,  "I -alpha". 


What  Confidence  Interval  do  you  want? 
(. 90,. 95, Your  Choice)?  .90 


We  had  decided  after  talking  to  the  boss,  that  an  alpha  of  ".10"  was  appropriate,  so  vve 
input  ".90  "  and  strike  the  "enter"  key.  After  the  machine  has  accepted  the  first 
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parameter,  it  will  return  with  the  following  screen  display  requesting  that  you  input  the 
desired  "population  proportion". 


What  is  the  population  proportion? 
(Q  Value)?  .95 


We  input  as  the  proportion  of  the  population  the  value  stated  in  the  contract 
specifications.  His  contract  states  that  95%  of  his  straps  will  meet  or  exceed  the  rated 
load  capacity,  therefore  we  input  ".95". 

The  machine  will  now  require  information  from  you  concerning  the  values  for  "R", 
and  ”M".  If  either  one  of  these  values  is  equal  to  zero,  then  the  tolerance  limit  is  one 
sided.  (R  =  0,  means  the  limit  is  only  upper  bounded,  while  an  M  =  0  is  lower  bounded.) 
The  most  commonly  used  values  for  R  and  M  are  1,  and  that  is  what  we  will  use  in 
this  example.  (  Note  -  The  sum  of  R  and  M  must  be  at  least  equal  to  "/"  or  you  will 
remain  in  the  loop.) 


While  the  inverse  chi  square  distribution  is  being  calculated,  the  following  self 
explanatory  screens  are  displayed. 
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PM'lflUHJ-IJVUM 


Upon  completi.  n  of  the  calculations,  the  screen  will  return  with 
following  output  data, 


To  be  .90  confident  you  need 
an  N  >  or  =  to  76. 3671875 


TO  CONTINUE  -  STRIKE  ANY  KEY 


From  this  outcome,  we  can  inform  our  boss  that  in  order  to  be  90%  confidence 
that  95%  of  the  contractors'  straps  are  indeed  exceeding  the  load  specification,  in  a 
given  lot. 

You  have  now  completed  the  Tolerance  limit.  We  strike  a  key  and  the  following 
screen  appears: 


Strike  'C1  to  do  another  test 
If  you  are  finished  strike  any  other  key 


If  you  strike  "C",  the  machine  will  reinitialize  itself.  Striking  any  other  key  followed  by 
the  "f8"  function  key  returns  you  to  the  M100  main  menu.  Do  not  strike  the  "enter" 
key  after  making  your  selection. 


|b.  |  Cochran  Test 

General 

The  Cochran  test  for  related  observations  is  an  attempt  to  apply  the  power  of  the 
analysis  of  variance  design  to  nominal  level  data.  The  test  layout  is  structured  with  the 
treatments  across  the  columns,  and  the  persons  or  blocks  down  the  rows.  The 
treatment  outcomes  are  dichotomous,  with  outcomes  only  indicated  as  "  1"  or  "0". 
Each  block  receives  every  treatment,  so  this  is  a  randomized  complete  block  design.  If 
there  are  only  two  treatments,  this  test  reduces  to  the  McNemar  test. 

*  *  Example  Problem  *  * 

Problem  Setting 

There  are  three  new  medications  for  sea  sickness  which  you  have  been  asked  to 
test  in  order  to  determine  which,  if  any,  should  be  stocked  on  board  Navy  vessels.  On 
I  one  aircraft  carrier  you  randomly  select  15  sailors  out  of  all  on  board  who  have  a 

|  history  of  sea  sickness.  You  further  subdivide  these  15  into  3  subgroups,  and  decide 

'  that  each  time  there  is  a  storm,  one  of  the  three  groups  will  receive  one  of  the  drugs, 
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until  each  group  has  been  administered  all  three  drugs.  This  design  is  to  overcome  the 
influence  of  storm  intensity.  If  drugs  work,  and  the  sea  sickness  abates,  then  they 
receive  a  T  otherwise  they  received  a  "0".  After  the  first  three  storms,  you  have  the 
following  data: 


SAILOR 


MEDICINE  USED 
B 


Total s 


Totals  8  11  8 

Let's  determine  if  there  is  any  statistically  significant  differences  between  the  three 
medications. 

Computer  Solution  to  Problem 

We  collected  nominal  level,  dichotomized  treatment  data  on  our  randomly 
selected  blocks,  the  sailors;  therefore,  we  will  select  "3"  from  the  Module  Main  Menu. 
(Sole  ihal  it  is  not  necessary  to  strike  the  "enter"  key  after  selecting  your  test.) 


1-Tolerance  Limits 
3-Cochran  Test 
5-Chi  Square  Test 


2-McNemar  Test 
4-Median  Test 
6-Goodness  of  fit 


PLEASE  MAKE  YOUR  SELECTION  <l-6>  3 


The  following  screen  will  now  appear: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  1 Y '  for  yes  or  1 N 1  for  no 


We  answer  with  a  Capital  "Y"  to  the  request  for  information  on  the  test  assumptions. 
A  Capital  "N"  would  have  caused  the  machine  to  skip  the  data  assumption  loop.  The 
following  screen  will  appear: 


You  selected  the  Cochran  Test 
this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


We  strike  any  key  to  continue. 

The  following  screen  now  will  appear: 


1- Randomly  selected  treatment  blocks 

2- Dichotomous  treatment  outcomes  (0,1) 

TO  CONTINUE  -  STRIKE  ANY  KEY 
If  you  fail  the  assumptions  then 
to  return  to  the  menu  input  'M 


The  data  we  have  collected  seems  to  meet  the  assumptions  of  the  test  so  we  "strike  any 
key"  and  begin  to  input  the  data  necessary  for  the  calculations.  If  you  make  an  error  in 
the  selection  of  the  test,  or  fail  to  meet  the  assumptions,  input  an  "M"  to  return  to  the 
Module  Main  Menu. 


First,  the  program  requests  information  on  the  number  of  columns  or  treatments. 
In  our  example,  three  different  drugs  were  administered,  so  we  input  a  "3". 


After  striking  the  "enter"  key,  the  M100  begins  a  loop  requesting  that  you  input  the 
column  totals.  The  loop  will  execute  three  times  for  our  example. 


We  input  8,  as  the  column  "1"  total,  11  as  the  column  "2"  total,  and  finally  8,  as 
the  column  "3"  total. 

The  machine  now  asks  about  the  number  of  "blocks"  that  we  have  in  the  design. 
We  had  15  sailors,  which  wc  randomly  selected  from  the  population  of  sailors  who 
suffer  from  sea  sickness.  {Note  -  there  are  no  size  limitations  on  the  number  of  rows  or 
the  number  of  columns  that  you  may  have  for  this  test,  except  those  due  to  the  MI00 
memory  size.) 


The  M100  now  enters  the  input  loop  for  the  row  totals.  For  our  example  data  we 
input  Fifteen  different  row  totals,  one  by  one,  as  the  Following  prompt  appears  15  times 
successively. 


Input  the  row  1  total?  2 


When  the  last  row  total  is  input,  the  M100  will  calculate  the  Cochran  test  statistic. 
The  following  self  explanatory  screen  will  appear  while  the  calculations  are  in  progress. 


When  completed  the  following  screen  will  appear: 


Your  test  statistic  is 
equal  to  1. 63636 
your  'Alpha  Hat1  is  .441233 

TO  CONTINUE  -  STRIKE  ANY  KEY 


Our  data  produced  a  "rest  statistic"  of  1.63636,  with  two  degrees  of  freedom  (C-l). 
The  small  number  of  degress  of  freedom  cause  "Alpha  Hat"  to  be  extremely  large.  We 
therefore  cannot  say  that  there  is  any  difference,  statistically,  between  the  three 
medications  we  tested.  We  might  recommend  to  the  boss  to  try  other  medicines,  to 
stock  the  cheapest  one,  or  take  "B"  since  numerically  it  did  perform  better. 


*,  N  %  %  -.  >  ■„ 


This  concludes  the  Cochran  test.  We  strike  a  key  and  the  following  screen 
appears: 


If  you  wish  to  perform  another  test,  strike  "C"  and  the  machine  will  reinitialize  itself.  If 
you  do  not  wish  to  continue,  then  strike  another  key,  and  followed  by  the  "f8"  key  to 
return  to  the  M100  main  menu.  You  need  not  strike  the  "enter"  key  after  making  your 
selection. 

|c.  |  Median  Test 

General 

The  median  test  is  used  to  test  whether  "n"  independent  samples  have  different 
medians,  or  central  tendencies  indicators.  It  is  really  a  special  application  of  the  chi 
square  test  with  the  row  marginals  fixed  by  the  value  of  the  "grand'  median.  (Note  - 
You  can  extend  this  test  to  examine  other  quantiles  by  simply  changing  the  criteria  for 
the  rows.)  [Ref.  2:p.  174). 

*  *  Example  Problem  *  * 

Problem  Setting 

Let  us  assume  that  the  Army  is  procuring  155mm  howitzer  shells  from  three 
different  manufacturers.  It  has  been  observed  when  units  are  firing  for  qualification 
that  their  CEP  (Circular  Error  Probable)  is  not  systematic.  You  have  been  assigned  to 
investigate  this  problem.  You  believe  that  it  is  possible  that  the  shells  from  the 
manufacturers  are  are  not  all  within  the  same  degree  of  error  tolerance.  You  collect 
data  on  the  next  1000  rounds  (ired  at  the  range,  directing  that  the  units  be  issued 
rounds  from  only  one  manufacturer  on  each  of  the  three  days  they  do  qualification 
firing.  In  this  way  you  will  attempt  to  control  for  the  possiblity  the  CEP  error  is  due  to 
unit  performance. 


fll 


Using  all  1000  rounds,  you  determine  the  "grand"  median  CEP  is  95  meters.  You 
break  out  the  projectiles  by  manufacturer,  and  the  data  looks  as  follows: 


MANUFACTURER 

Totals 

A 

B 

C 

>95 

287 

105 

108 

500 

<=95 

47 

227 

226 

500 

Total s 

334 

332 

334 

Computer  Solution  to  Problem 


The  measurement  of  CEP  errors  in  the  field  is  not  extremely  precise,  so  we  will 
only  consider  this  data  to  be  ordinal.  We  are  interested  in  testing  if  the  CEP  for  the 
different  manufacturers  projectiles  are  the  same.  From  the  Module  Main  Menu  we 
select  "4".  ( Note  that  it  is  not  necessary  to  strike  the  "enter"  key  after  selecting  your 
test.) 

mInu 

1-Tolerance  Limits  2-McNemar  Test 
3-Cochran  Test  4-Median  Test 

5-Chi  Square  Test  6-Goodness  of  fit 

PLEASE  MAKE  YOUR  SELECTION  <l-6>  4 

The  following  screen  will  now  appear: 


We  will  check  our  assumptions  in  this  example  problem,  so  we  input  a  to  the 
request  for  information  on  the  test  assumptions.  The  next  two  screens  will  appear: 
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We  believe  that  we  have  data  which  meets  the  above  assumptions;  therefore,  we  strike 
any  key  and  continue.  If  you  discover  that  you  fail  to  meet  the  assumptions,  input  an 
"M". 


In  order  to  perform  the  calculations,  the  Ml 00  requires  that  you  input  the  number 
of  rows  and  columns  in  the  contingency  table.  We  will  input  a  "2"  for  the  number  of 
rows  and  "3“  for  the  number  of  columns,  since  we  had  three  munitions  manufacturers, 
striking  the  "enter"  key  after  both  inputs.  ( Note  -  the  median  test  will  only  accept  "2" 
rows  as  an  input,  and  these  must  be  input  so  that  the  "above  the  median  scores"  are  placed 
in  the  first  row.) 


Please  enter  the  number  of  rows  in 
the  contingency  table?  2 
Now  enter  the  number  of  columns?  3 


The  Ml  >0  now  begins  a  loop  requesting  that  you  input  the  cell  frequencies  begining 
with  the  first  row  first  column,  and  working  across  the  columns. 


Input  the  number  in  the  cell  in 

ROW-  1 
COLUMN-  1 


THE  NUMBER  ?  287 


We  input  the  numbers  287,  105,  108,  for  the  first  row,  and  47,  227,  226  for  the 
second  row.  In  the  example,  we  selected  the  2  row  totals  were  exactly  equal.  If  they 
had  not  been  equal,  the  M100  would  have  queried  you  to  determine  if  you  have  made 
an  error.  If  you  have  not,  simply  input  "N"  and  continue,  otherwise  you  have  anc'her 
opportunity  to  input  the  data.  After  the  last  number  is  entered,  the  following  screen 
will  appear. 


When  the  calculations  are  finished,  the  M100  will  return  the  following  screen 
displaying  the  test  results. 


Your  test  statistic  is 
equal  to  258.975 
your  'Alpha  Hat1  is  0 

TO  CONTINUE  -  STRIKE  ANY  KEY 


Here  our  observed  "test  statistic"  was  so  large  that  the  "Alpha  Hat"  was  smaller 
than  .00001,  and  has  been  rounded  to  zero.  We  can  say  that  there  is  statistical 
difference  between  the  rounds  manufactured  by  the  three  different  companies. 


This  concludes  the  Median  test.  We  strike  a  key  and  the  following  screen 


appears: 


Strike  'C'  to  do  another  test 
If  you  are  finished  strike  any  other  key 


The  machine  will  reinitialize  itself  if  you  strike  "C .  If  you  do  not  wish  to  perform 
another  test  then  strike  any  key,  and  then  the  "f8"  key.  This  returns  you  to  the  M100 
main  menu.  You  need  not  strike  the  "enter"  key  after  making  your  selection. 

|d.  |  Chi  Square  Test  for  Differences  in  Probabilities 

General 

The  Chi  Square  test  is  used  to  determine  if  the  the  probability  of  discrete 
categories  of  occurrences  are  the  same  in  several  populations.  It  is  also  used  is  to  test 
if  several  independent  samples  can  be  assumed  to  have  come  from  the  same 
population.  Although  this  test  is  used  very  frequently,  it  is  not  recommended  in 
situations  where  the  data  is  clearly  interval  level,  and  the  assumptions  of  the  F  test  are 
satisfied. 

*  *  Example  Problem  *  * 

Problem  Setting 

Suppose  that  you  work  for  the  Army  recruiting  command  at  a  recruiting  station 
located  in  a  metropolitan  area.  There  are  four  high  schools  located  in  the  community: 
a  private  high  school,  two  public  high  schools,  and  a  technical  trade  school.  In  the 
past,  you  have  had  the  recruiters  present  the  same  speech  at  all  the  high  schools  in  an 
attempt  to  convince  the  students  to  enlist  upon  graduation.  The  effects  of  the  "dog 
and  pony  show"  on  enlistments  docs  not  seem  to  have  been  the  same  in  all  the  high 
schools.  You  believe  that  perhaps  the  speech  should  be  varied  depending  on  the  type 
of  recruits  you  can  expect  to  get  from  each  school.  You  have  records  on  all  the 
enlistments  from  this  recruiting  station  for  the  last  two  years  and  these  records  include 
the  high  school  from  which  the  enlistee  graduated.  You  would  like  to  determine  if  the 
enlistment  profiles  by  MOS  (military  occupation  status)  vary  by  high  school.  You 
decide  to  divide  the  jobs  into  three  major  areas,  combat  arms  (infantry,  armor,  etc.), 
combat  support  arms  (engineers,  maintenance,  military  police,  etc.),  and  service 
support  arms  (finance,  legal,  medical,  etc.)  and  to  take  a  random  sample  of  all  the 


local  enlistments.  Your  null  hypothesis  is  that  there  is  no  difference  between  the  high 
school  enlistment  profiles. 


SCHOOL  TYPE 

Combat 

MILITARY 
Combat  Support 

SPECIALITIES 

Service  Support 

Private 

7 

14 

22 

Public  (1) 

18 

9 

6 

Public  (2) 

21 

11 

8 

Technical 

8 

23 

11 

Let's  solve  this  problem  using  the  M100. 

Computer  Solution  to  Problem 


We  have  frequency  count  data  from  three  independent  samples,  therefore  we  are 
limited  (assuming  we  do  not  violate  any  test  assumptions)  to  nominal  tests  for  "k" 
independent  samples.  Since  we  do  not  have  a  treatment/block  design  needed  to  perform 
a  Cochran  test,  it  would  seem  that  the  most  appropriate  test  on  this  module  is  the  Chi 
Square  test.  We  select  "5"  from  the  Module  Main  Menu.  ( Note  that  it  is  not  necessary 
to  strike  the  "enter"  key  after  selecting  your  test). 


MENU 

1-Tolerance  Limits  2-McNemar  Test 
3-Cochran  Test  4-Median  Test 

5-Chi  Square  Test  6-Goodness  of  fit 

PLEASE  MAKE  YOUR  SELECTION  <l-6>  5 


The  following  screen  will  now  appear: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  'Y1  for  yes  or  'N1  for  no 


\ 


We  believe  that  we  meet  all  the  assumptions,  but  we  will  check  the  assumptions  to 
ensure  that  we  have  not  overlooked  any.  We  input  a  Capital  "Y"  to  the  request  for  this 
information.  If  you  want  to  skip  the  loop  for  the  assumptions,  then  a  Capital  "N" 
would  have  been  the  appropriate  response.  Since  we  input  a  "  T'  the  following  two 
screens  will  appear: 


You  selected  the  Chi  Square  test 

this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


1- Independent  random  samples 

2- Each  observation  is  classified  into 
only  one  row  and  one  column 

TO  CONTINUE  -  STRIKE  ANY  KEY 
If  you  fail  the  assumptions  then 
to  return  to  the  menu  input  M1 


We  know  that  we  have  random  samples,  and  that  a  recruit  could  have  graduated  from 
only  one  of  the  high  schools  in  question.  We  strike  key  and  continue.  If  you  discover 
that  you  fail  to  meet  the  assumptions,  input  an  " \t" .  This  input  returns  you  to  the 
Module  Main  Menu. 

The  calculations  require  that  you  input  the  number  of  rows  and  columns  in  the 
contingency  table.  We  will  input  a  "4"  for  the  number  of  high  schools  shown  in  the 
rows  and  "3"  for  the  military  specialties  shown  in  the  columns,  striking  the  enter  key 
after  each  input. 


Please  enter  the  number  of  rows  in 
the  contingency  table7  4 
Now  erter  toe  number  of  columns?  3 
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The  M100  now  begins  a  loop  requesting  that  you  input  the  cell  frequencies 
begining  with  the  first  row,  first  column  and  working  across  the  columns. 


Input  the  number  In  the  cell  in 

ROW-  1 
COLUMN-  1 

THE  NUMBER  ?  7 


We  input  all  the  data  as  it  is  laid  out  in  the  "problem  setting"  portion  of  this 
discussion.  When  the  last  row  and  column  data  has  been  input,  the  following  screen 
will  appear: 


The  machine  will  calculate  the  test  statistic  from  the  data,  and  when  it  has  finished 
it  will  return  the  following  screen  displaying  the  results  for  the  problem  data  input. 


Your  test  statistic  is 
equal  to  30.0849 
your  'Alpha  Hat1  is  3. 8E-05 

TO  CONTINUE  -  STRIKE  ANY  KEY 


With  these  results,  wc  should  reject  our  null  hypothesis.  It  would  appear  in  this 
ease,  that  there  does  seem  to  be  a  difference  in  the  enlistment  profiles  of  these  high 
schools.  I  he  recruiters  might  find  their  time  more  wisely  invested  if  they  spoke  about 
v-omhat  arms  opportunities  at  some  schools,  and  combat  support  arms  or  service 
'upp  rt  arms  opportunities  at  others. 


This  concludes  the  Chi  Square  test.  We  strike  a  key  and  the  following  screen 
appears: 


If  you  wish  to  perform  another  test  strike  "C'  and  the  machine  will  reinitialize  itself. 
To  return  to  the  M100  main  menu  strike  any  key,  and  followed  by  the  "J8“  key.  You 
need  not  strike  the  "enter"  key  after  making  your  selection. 

|e.  IMcNemar  Test 

General 

The  McNemar  test  evaluates  independent,  bivariate,  nominal  level,  random 
variables  from  related  samples  in  order  to  determine  if  we  can  ".  .  .detect  a  difference 
between  the  probability  of  (0,1)  and  probability  of  (1,0)"  [Ref.  2:p.  130).  Most 
frequently  this  test  is  used  to  determine  the  significance  of  an  intervening  variable.  For 
example,  when  we  take  a  before  and  after  sample,  we  are  attempting  to  determine  if 
some  event  between  the  two  sample  measurements  had  any  significant  impact.  (Note  - 
even  though  we  obtain  significance  we  must  attempt  to  ascertain  if  the  intervening 
variable  measured  is  really  the  correct  one.) 

*  *  Example  Problem  *  * 

Problem  Setting 

Assume  that  you  are  the  training  officer  of  a  basic  training  unit.  You  have 
received  some  films  which  are  reported  to  be  effective  in  changing  attitudes  about 
drugs.  You  view  these  films,  and  feel  that  the  drug  users  in  the  films  are  slightly 
glamorized.  You  believe  that  this  might  make  these  films  less  effective  on  moderating 
drug  use.  You  decide  to  test  your  belief  by  taking  a  group  of  viewers  "before"  the  film 
and  asking  them  questions  about  drug  use.  After  the  film  series,  you  ask  the  same 
questions  to  all  .hose  in  the  original  group  who  viewed  the  entire  scries.  The  results  on 
one  such  question  arc  as  follows: 


Would  you  use  marihuana? 


(AFTER  FILMS) 

(BEFORE  FILMS) 

Yes 

No 

Yes 

18 

12 

No 

8 

15 

Let's  test  our  original  premise  by  conjecturing  that  there  will  be  no  significant 
difference  in  the  (0,1)  and  (1,0)  shifts. 


Computer  Solution  to  Problem 


We  have  related  samples  with  a  dichotomous  response  to  a  question.  The  level 
of  our  data  is  nominal,  therefore  it  would  seem  that  the  most  appropriate  test  on  this 
module  is  the  McNemar  test.  We  select  "2"  from  the  Module  Mam  Menu.  (Note  that 
it  is  not  necessary  to  strike  the  "enter"  key  after  selecting  your  test.) 


MENU 

1-Tolerance  Limits  2-McNemar  Test 
3-Cochran  Test  4-Median  Test 

5-Chi  Square  Test  6-Goodness  of  fit 

PLEASE  MAKE  YOUR  SELECTION  <l-6>  2 


The  following  screen  will  now  appear: 


Do  you  wish  to  check,  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  1 Y 1  for  yes  or  ‘N1  for  no 


We  will  check  the  assumptions.  We  input  a  Capital  "Y"  to  the  request  for  this 
information.  If  you  want  to  skip  the  loop  for  the  assumptions,  then  a  Capital  " S" 
would  have  been  the  appropriate  response.  Since  we  input  a  "Y'\  the  following  three 
screens  will  appear: 
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( Note  -  the  McNemar  test  which  uses  a  smaller  "n"  is  not  calculated  here.) 


s 

& 


If  you  fail  the  assumptions  then 
to  return  to  the  menu  input  M' 


We  know  that  the  set  of  "X"  responses  (before)  and  "Y"  responses  (after)  are 
independent  nominal  level  responses,  and  that  "N"  is  greater  than  twenty.  We  have, 
therefore,  met  the  assumptions  and  strike  any  key  and  continue.  If  you  discover  that 
you  fail  to  meet  the  assumptions,  input  an  ".V/".  This  input  returns  you  to  the  Module 
Main  Menu. 


We  input  the  four  cell  frequencies  shown  in  the  "problem  setting"  portion  of  this 
discussion.  When  the  last  row  and  column  data  has  been  input,  the  following  screen 
will  appear: 


The  machine  will  calculate  the  test  statistic  from  the  data,  and  when  it  has  finished 
it  will  return  the  following  screen  displaying  the  results  for  the  problem  data  input. 


From  the  results  of  the  data  you  should  retain  your  null  hypothesis.  However, 
even  more  importantly  you  should  inform  your  boss  that  this  film  has  an  adverse  effect 
on  attitudes  about  drugs  and  should  probably  be  discontinued. 

This  concludes  the  McNcmar  test.  We  strike  a  key  and  the  following  screen 
appears: 


If'  you  wish  to  perform  another  test  strike  C"  and  the  machine  will  reinitialize  itself 
lo  return  to  the  Mint)  man.  menu  strike  am  key,  and  followed  b\  the  ' j\  key  You 
need  not  strike  the  enter  key  alter  making  your  selection. 
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|f.  |  Goodness  of  Fit  test 

General 

The  Chi  Square  Goodness  of  Fit  test  allows  a  single  sample  to  be  compared  with  a 
theoretical  distribution  to  determine  whether  the  sample  could  have  been  drawn  from 
this  hypothesized  population. 

*  *  Example  Problem  *  * 

Problem  Setting 

Assume  that  as  an  analyst,  you  have  conducted  a  survey.  You  are  not  certain  of 
the  underlying  population  distribution  and  want  to  use  only  nonparametric  methods  to 
analyze  all  the  data.  Your  boss,  however,  would  like  you  to  use  parametric  tests.  You 
decide  to  check  your  data  for  normally,  since  this  is  an  underlying  assumption  for 
many  parametric  tests.  You  have  the  following  results  on  the  ages  of  your 

respondents: 

17  18  18  18  18  19  19  19  19  20  20  20  20  21  23 

24  25  26  28  29  30  30  34  37  39  40  42  47  48  50 

We  calculate  the  mean  and  the  standard  deviation  for  this  set  of  ages.  Our  mean  age  is 
27.3  years  and  the  standard  deviation  is  10.2  years.  Let's  check  to  see  if  this  set  of  data 
are  consistent  with  a  normal  population. 

Computer  Solution  to  Problem 

We  wish  to  do  a  goodness  of  fit  test,  so  we  select  "6"  from  the  Module  Main 
Menu.  (Sole  that  it  is  not  necessary  to  strike  the  "enter"  key  after  selecting  your  test.) 

MENU 

1-Tolerance  Limits  2-McNemar  Test 

3-Cochran  Test  4-Median  Test 

5-Chi  Square  Test  6-Goodness  of  fit 

1 

PLEASE  MAKE  YOUR  SELECTION  <l-6>  6  | 


The  following  screen  will  now  appear: 


Do  you  wish  to  check  the  assumptions 
of  the  test  that  you  have  selected? 

Just  input  1 Y 1  for  yes  or  ' N 1  for  no 


We  will  check  the  assumptions.  We  input  a  Capital  "Y"  to  the  request  for  this 
information.  If  you  want  to  skip  the  loop  for  the  assumptions,  then  a  Capital  "N" 
would  have  been  the  appropriate  response.  Since  we  input  a  “Y‘  two  screens  will 
appear: 


You  have  selected  the  Goodness  of  fit 
this  test  assumes  you  meet  the 
following  assumptions 

TO  CONTINUE  -  STRIKE  ANY  KEY 


1- An  independent  random  sample 

2- Data  level  appropriate  to 

the  distribution  selected 


TO  CONTINUE  -  STRIKE  ANY  KEY 
If  you  fail  the  assumptions  then 
to  return  to  the  menu  input  ' M 1 


We  drew  a  random  sample  when  we  administered  our  questionnaire,  and  our  data  on 
age  is  adequate  to  permit  the  calculation  of  the  two  parameters  necessary  for  a  normal 
distribution.  We  therefore  believe  that  we  meet  the  assumptions  of  the  test.  We  strike 
any  key  and  continue.  If  you  discover  that  you  fail  to  meet  the  assumptions,  input  an 
".V/".  This  input  returns  you  to  the  Module  Main  Menu. 


The  M100  now  requests  that  you  indicate  the  theoretical  distribution  against 
which  you  wish  to  test  your  sample  data. 
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Select  your  theoretical  distribution 

1-  Normal  Distribution 

2-  Uniform  Distribution 

3-  Poisson  Distribution 

PLEASE  MAKE  YOUR  SELECTION  <l-3>  1 


We  input  a  "\"  as  our  choice,  since  we  assumed  a  normal  distribution.  (Note  -  do 
not  strike  the  "enter"  key  after  making  your  selection.) 

The  M 100  will  now  request  information  on  the  number  of  discrete  categories  into 
which  you  wish  to  break  the  data.  (Note  -  this  would  also  be  required  if  you  were 
performing  a  Poisson  or  a  uniform  goodness  of  fit  test.  The  classes  on  the  Poisson 
distribution  are  assumed  to  run  from  0  to  one  less  than  the  number  of  classes,  because 
the  Poisson  runs  from  0  to  infinity.  Therefore  the  M 100  adjusts  the  number  of  classes 
you  input  to  insure  that  this  occurs.  For  this  example,  we  chose  to  "discretize"  the 
normal  curve  into  6  sections.  How  many  are  chosen  is  a  matter  of  judgement  on  the 
part  of  the  analyst. 


Please  input  the  number  of  classes?  6 


We  strike  the  "enter"  key.  The  next  screen  appears  requesting  first  a  mean,  and 
then,  without  the  screen  clearing,  a  variance.  (For  the  Poisson  distribution,  the  M100 
requests  that  lambda  be  input,  while  for  the  uniform,  this  screen  is  skipped.)  We  may 
estimate  these  parameters  from  our  data,  or  we  may  use  known  values  for  both  or 
either.  The  effect  of  calculating  both  parameters  from  the  sample  data  will  be  a  loss  of 
two  degrees  of  freedom.  In  our  example  we  estimated  the  parameters  from  the  sample 
data. 
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Based  on  the  sample  mean  and  the  variance  the  program  calculates  the  scores  t-o 
be  used  as  the  class  boundaries.  (When  you  test  for  a  Poisson  or  a  uniform,  this  step 
simply  requests  the  input  of  the  number  of  cases  in  each  classes.) 


Input  the  #  of  observed  data  points 
<  or  =  17.3357  ?  1 


Input  the  #  of  observed  data  points  that 
are  >  17. 3357  and  <  or  =  22. 8678  ?  13 


The  number  of  scores  falling  in  each  of  the  classes  are  input.  In  our  example,  we 
input  successively  1,  13,  4,  4,  2,  and  6.  When  we  strike  the  "enter"  key  after  the  last 
input  the  following  screen  appears  requesting  information  of  the  parameters. 


N'! 

How  many  parameters  did  you  estimate 

' 

N.’ 

(0,  1  OR  2)?  2 

When  the  Poisson  distribution  is  tested  however,  the  M100  will  ask  if  you  estimated 
lambda.  The  estimation  of  this  parameter  will  cause  the  loss  of  one  degree  of  freedom. 

The  next  self  explanatory  screen  appears: 


V 


The  machine  will  calculate  the  test  statistic  from  the  data,  and  when  it  has  finished 
it  will  return  the  following  screen  displaying  the  results  for  the  problem  data  input. 

Your  test  statistic  is 
equal  to,  18.4 
your  'Alpha  Hat*  is  3. 64E-04 

TO  CONTINUE  -  STRIKE  ANY  KEY 

From  these  results,  you  cannot  perform  parametric  tests  on  the  data  which 
require  the  underlying  assumption  of  normality. 

This  concludes  the  goodness  of  fit  test.  We  strike  a  key  and  the  following  screen 
appears: 


If  you  wish  to  perform  another  test,  strike  "C"  and  the  machine  will  reinitialize  itself. 
To  return  to  the  M100  main  menu  strike  any  key,  and  followed  by  the  "/T  key.  You 
need  not  strike  the  "enter"  key  after  making  your  selection. 


APPENDIX  C 

MODIFYING  PROGRAM  MODULES 

To  redimension  the  arrays  for  the  BINOM  module  change  the  dimensions  of  all 
the  arrays  listed  in  line  130,  shown  below. 

130  DIM  CDF!  (30), D!  (30)  .PDF!  (30) 

The  line  to  be  changed  in  the  CHISQ  module  is  line  20  shown  below. 

20  DIM  X! ( 15, 15) ,C! (15), R! ( 15) ,N(3) ,N! (3) 

In  this  line  you  need  not  change  all  the  arrays  however,  only  the  first  three  variables 
X!,  C!,  and  R'. 


APPENDIX  D 

POSSIBLE  PROGRAMMING  PROBLEMS 


My  keyboard  entries  are  not  accepted. 

This  usually  is  solved  by  checking  to  see  that  the  "Caps  Lock"  key  is  depressed.  If 
this  does  not  help  the  problem  you  are  attempting  to  input  a  parameter  which  is  not 
within  the  bounds  allowed.  Check  the  sample  problem  in  the  user's  manual. 

When  I  try  to  run  the  program  I  get  an  "?OM  Error". 

This  error  will  occur  if  you  do  not  have  adequate  memory  available  to  initialize  the 
arrays,  f  irst,  check  to  see  that  the  "paste"  buffer  does  not  contain  any  information.  If 
it  is  already  empty,  then  it  will  be  necessary  to  erase  some  programs.  To  do  so,  enter 
basic  and  type  the  following: 

KILL  "<program  designatfor?>"  <enter> 

1  ran  the  program  and  it  returned  a  "7BS  error". 

You  will  receive  a  bad  subscript  error  if  you  exceed  the  dimensioned  size  of  the 
arras  l!  this  occurs,  then  see  the  appendix  on  redimensioning  arrays,  or  solve  the 

problem  by  hand. 

I  he  M  100  skipped  a  screen. 

I  here  are  two  possible  causes  of  this,  f  irst,  you  struck  the  keyboard  too  hard 
a  her;  putting  a  parameter.  Second,  you  struck  the  <  enter  >  key  when  it  was  not 

I  tried  to  copy  the  program  into  BASK  and  got  a  "?DS  error". 

I  :;;s  i.si.alh  indicates  that  there  are  carriage  returns  imbedded  in  program  i.nes 
1  Lew  ma\  ue  removed  K  scrolling  through  the  program  to  find  them,  or  using  the 
p:  c:  tm  <  R<  H  1  1\1  lrom  the  M  l"<>  disk  to  remove  them. 
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APPENDIX  E 

BINOMIAL  DISTRIBUTION  PROGRAM 


100  CLS 

120  QEFSNGA: DEFINTM: OEFINTN;  QEFINTK:  DEFINTT: 

130  DIM  CDF!  (30) ,  D!  (  30) ,  PDF!  (  30 ) 

140  $=1 

160  GOSU6  10020  'THIS  IS  BINOMIAL  MENU  SUBROUTINE 
ISO  CIS:  PRINT;  PRINTTABf  3)  ''Do  you  wish  to  check  the 
135  PRINT:  PRINT;  PRINTTAB(3),  "of  the  test  that  you  h 
190  PRINT@203,  "Just  Input  'Y'  for  yes  or  rNT  «« 
200  e$=lNKEY5:IF  B$=  "  "GOTO  200 
220  IF  B$-  Y"  OR  8$="y"  GOTO  280 
240  IF  3$="N"  OR  5§="n"  GOTO  1000 


OE  F I  NT  X 


190  PRINTO203 ,  "Just  Input  V  for  yes  or  'N 

200  e$=lNKEY5:IF  B$=  "  "GOTO  200 

220  IF  B$-  Y"  OR  8$="y"  GOTO  280 

240  IF  3$="N"  OR  BS="n"  GOTO  1000 

260  CLS:  PRINT(3125 ,  The  response  must  be  '  Y '  or 

270  PRINT(3174  "TRY  AGAIN":  GOTO  180 


for  no 


assumptions" 
^ve  selected’" 


280  GOSUB  110(50 

1CC0  ON  QT  GOTO  2000.2100.2700,3200,2700 

2000  BINOMIAL  TE$T 

2020  GOSUB  18000  YSUBROU?lNE  TQ  INPUT  AL 


2040  GOSUB  15800  '  PDF  CALC 

2C60  GOSUB  16080  '  ALPHA 

2065  GOSUB  13020  '  INKEY 

2070  GOTO  13200  'ALPHA  HATS 

2100  Y****  QUANTILE  TEST 


PDF  CALCULATIONS 

ALPHA 

INKEY 


ALL  PARAMETERS 


2111  GOSUB  18000  rlNPUT  VARIABLES 
2121  GOSUB  15800  'PDF  CALCULATION! 


2121  GOSUB  15800 
2131  GOSUB  16080 
2141  GOSUB  13740 


'PDF  CALCULATION! 
'ALPHA 

'OUTCOME  OUTPUTS 


2151  ON  T  GOTO  2161 .2161 ,2181 

2161  PRINTTAB(2);nttlf  is  less  than  or  equal  to  LN 

2171  IF  T=2  THEN  2201 

2181  IF  UN+1>N  THEN  UN=UN-1 

2191  PRINTTAB(2); "'T2'  is  greater  than  ";UN 

2201  GOSUB  13020 


2191  PRINTTAB(2); n'T2‘  is  greater  than  ";UN 
2201  GOSUB  13020 
2211  GOTO  13200 

2700  Y***“  SIGN  TEST/COX  STUART 

2720  '  PL=#  OF  (♦VS 

2740  '  MI=#  OF  (-)'S 

2780  ’  N=SUM  OF  (♦)  AND  (-)'S 

28J0  '  S  IS  A  LOGICAL  COUNTER  FOR  SIGN  TEST 

2820  CLS: PRINTP120,"  Input  the  #  of  times  X<Y  (♦)"; 

2340  INPUT  PL: CLS  ,  „ 

2860  PRINT0120 . "  Input  the  #  of  times  X->Y 
2380  INPUT  MI: OlS 
2900  P=.  5: S=1 
2920  N=PL+MI 

2940  GOSUB  15320  'INPUT  ALPHA  TOTAL 
2960  GOSUB  15480  'INPUT  TAIL  VALUES 
2980  S=T+1 

3000  GOSUB  15800  'CALCULATES  PDF, AND  CDF 

3020  GOSUB  16740  'ALPHA 

3025  GOSUB  13020  'INKEY 

3040  GOTO  13200  'ALPHA  HATS 

3200  Y****  QUANTILE/C  I 

3220  GOSUB  18000  'THE  INPUT  VARIABLES 

3240  GOSUB  15800  ’PDF  CALCULATIONS 

3260  G0SU8  16080  1  THE  ALPHA 

3300  CLS: PRINT: PRINT"  The  extreme  tail  values  are  presented 
3305  PRINT"be1ow.  You  will  need  to  select  the  N's” 

3010  PRINT"  which  will  determine  the  calculation' 

3015  PRINT"  of  your  total  ALPHA. " 

3320  PRINT: PRINTTAB(7)-/'TO  CONTINUE-STRIKE  ANY  KEY" 


3340  K$=INK£Y$: IF  K$  =  "" 
T  GOTO  3380,3660, 


40 


3340 


3360  ON 
3380  CL' 

3385  P«INT''08S#^TAB(5);"L0WER  COF";  TAB(21);  "OBS#";  TAB( 28);  "UPPER  1-CDF" 

J  ‘  "LSI 
LSI 


3390  FQ 
3400  IF 
3420 
3430 


00  IF  LN=0  \ 
20  IF  UN=N  T 
30  ZZ  =  LN  > 


Th||N  N2=fx4j  ^ 


.  .  rTAB(5k  ALPHA 

3540  if  “  :  r: 

3560  PR'NTTAb(4) 

3565  IF  UN-*- 1  '  H  . 

3580  IF  T=P  THEN  3620 
3600  PRINTTAB(4),  i1YT2 
3620  IF  T=1  THEN  3820 
3630  GOSUB  13020 
3635  GOTO  13100 


N$X-§) 

TAB( 21 ); 

Input  an  N  for  lower  OBS#"; LN 
t  an  N  for  upper  OBS#"; UN 
:  H2=1-CDF! (UN+1): 


3440  PRINT  ZZ;  TAB(5);  COF!  (ZZ+l);  TAB(21);  UN+N2;  TAB( 30 ) ;  1-(CDF!  (UN+N2+1)) 

3460  PRINT, INPUT" 

3480  INPUT1'  Inc 
3500  H1=CDF!(IN- 
3520  CLS:  PRInTTaL  - 
“  IF  r=3  r-OTO  ^660 

; "YT 1 1  lower  C.  I. 

GOTO  3640 


Is  equal  to";HT:  PRINT 
Is  <  or  =  ";  LN-*- 1 


upper  C.  I.  Is  >  or  =  ";UN+1 


3640  UN=IJN-1:  GOTO  3600 

3660  CLS: PRINTTAB(9); wOBS#": TAB(19): "LOWER  CDF" 

3680  FCR  X  =  1  TO  3:  IF  LN=0  T&EN  N1  =  (S<-1)  ELSE  Nl=(X-2) 

3700  PRINTTAB(9)  LN+N1; TAB( 19). CDF! ( LN+N1+1): NEXTX 

3720  PRINT: INPUT"  Input  the  N  tor  the  lower  OBS#"; LN  : HT=CDF! ( LN+1 ) 

3730  GOTO  3520 

3740  CLS  PRINTTAB(9); "OBS#", TAB( 19);  "UPPER  1-CDF" 

3760  FOR  X=1  TO  3:  IF  UN=N  T^EN  N2=U-3)  ELSE  N2=(X-2) 

3780  PR INTTABf 9 )  UN+N2, TAB( 19); 1-CDF? (UN+N2+1): NEXTX 
1300  PRINT: INPUT  "  Innut  the  N  for  the  upper  OBS#  ; UN 
3310  HT  =  l:CDF'JiJN^l):  GOTO  3520 


3320  GOSUB  13020 
3340  CLS:  PRINT3120 ," I f  you  want  the  data  value  at  T1  AND  T2": 

3845  PRINT0215 ." 1 nput  an  V" 

3350  PRINTTAB(3);  otherwise  any  other  letter  "; 

3360  INPUT  XV  IF  X$<>"Y"  GOTO  13100 
3830  GCSUB  13120 
39U0  GOS  ;B  13300 

3  320  OlS:  PRINT:  PR  I  NT  T  AB(  10);  "  T 1  =  Q!  (  LN+1 );  "  and  T2  =  ";D!(UN+1) 

PRcNT 

5)4)  GCSjB  13020 


:o 


uvT0 13100 

I  *  ft  *  * 


;020  PR INTTABf 18). "MENU 
022  PR  I  NT" 1-B i nomi al  rest" 
'025  PRJNT"  JrCox-S^uart  tes 


SUBROUTINF  MENU  FOR  THE  BINOMIALLY  DISTRIBUTED  TESTS 

’  a  \  »  m  c  w  m  » 


io : 


5-Sign  Test" 

PLEASE  MAKE  YOUR  SELECTION 
-  GOTO  10060 


:  TAB( 20 )- "2-Quant  1 le 
l"; TAB(20), "4-Quanti 1 


Test" 
e  C.  I. 

1-5'"; 


111 
1  1245  PR  IN 
11^50  -RIN 


i  -i 
1  1 


PRINTTAB(IO), 

PRIST:  PRINT."  . 

's-INKE^S’ IF  K$  =r 
;T=VALt<S) 

N  Qr-o  GOTO  10240 
IF  CM  GOTO  10240 
PJRN 

S:  PRINT0123 "vou  nave  made  an  improper  selection" 
PR INTTABf 1 3 ): Ato  try  aga’n" 

PR  INTTABf  7 )  'TO  CONTINUE  -  STRIKE  ANY  KEY" 
^S=INKE<S: *$="*  GOTO  10280 
'1 1_ S :  GOTO  10020 
REMN 

- -  SUBROUTINE  FOR  THE  ASSUMPTIONS  OF  THE  TESTS 

ON  G*  GOTO  11040,11260,11340,11260,11300 
S:  r?IN • 'j’34  , " You  selected  the  Binomial  Test" 
f6l" this  test  assumes  you  meet  the' 

J; "following  assumptions" 

B  13020 


,60  PRINT: 


'AB 

TAB 


m 


PRINT343 . " 1-On 1 y  Nominal  Level  Data"; 
11100  PRINT033 , "2-  P *  Constant  Over  All  Trials’ 


K4 


11105 

11120 

11140 

11160 

11180 

11200 

11220 

11240 

11260 

11265 

11270 

11280 

11290 

11300 

11305 

11310 

11320 

11325 

11330 

11335 

11340 

11345 

11350 

11360 

11370 

13000 

13020 

13040 

13060 

13100 

13110 

13120 

13140 

13160 

13175 

13200 

13210 

13215 

13220 

13225 

13230 

13240 

13260 

13280 

13300 

13320 

14000 

14020 

14040 

14080 

14100 

14110 

14120 

14140 

14160 

14130 

14190 

14195 

14200 

14220 

14240 

14260 

14300 

14320 

14340 

14360 

14380 

14400 

14420 


PRINTTAB(3)"3-A1 1  Observations  Are  Independent" 
PRINT@203.jT0  CONTINUE  -  STRIKE  ANY  KEY  * 

PRINT0243,' If  you  fall  the  assumptions  then" 
PRINT@283,"to  return  to  the  menu  input  rMT; 

E$=INKEY$: IF  E$='™  GOTO  11180 
IF  E$<>  M  GOTO  1000 
CLS: GOTO  160 
RETURN 

CLS:  PRINT(385."You  have  selected  the  Quantlle/CI  " 
PRINTTAB(6j: "thi s  test  assumes  you  meet  the" 

PRINT! 4BL.10V,  .following  assumptions":  PRINT:  GOSUB  13020 
CLS:  PRINT(38l  "1-OrdinaT  or  Higher  Level  Data" 

PRINTTABC 3); a2-A1 1  Observations  are  IID  ": GOTO  11120 
CLS: PRINT@82LhYou  have  selected  the  Sign  Test" 
PRINTTAB(5k  "this  test  assumes  you  meet  the" 
PRlNTTABf9V'fol lowing  assumptions":  PRINT:  GOSUB  13020 
CLS:  PRINT:  PRINTTAB( 3); "1-Orai nal  or  Higher  Level  Data" 
PRINTTAB(3);rt2-Paired  Observations1' 


PRINTTAB( 3); "2- 
PRINTTAB(3j!"3- 
PRINT0163.^-Ca 
CLS: PRINT$85  ."Y 


Outcomes 


Classes 


ases  IID  Within  Classes 


f+,0R 

. 6oto  ii; 


veduvAlue 


PRINTTAB(6): "thi s  test  assumes  you  meet  the" 
PRINTTAB(.10V,  "following  assumptions":  PRINT:  GOSUB  13020 
CLS:  PR  I N  T(383 .  1-Ordinal  or  Higher  Level  Data" 
PRINTTAB(3); rt2-All  Observations  Are  IID": GOTO  11120 
*****  SUBROUTINE  INKEYS 

PRINTTABI7}- "TO  CONTINUE  -  STRIKE  ANY  KEY"; 

K$=INKEY$:  K$=it"  GOTO  13040 

RETURN 

CLS: PRINT@83 .  "Strike  'C'  to  do  another  test" 

PRINTP120 ."If  you  are  finished  strike  any  other  key" 
K$=INKEY$:  IF  K$=""GOTQ13120 
IF  K$-  C"  GOTO  13175 
CLS:  END 
CLS: GOTO  160 

'  ****  REQUEST  FOR  ALPHA  HATS 
*  VARIABLES  AT=ALPHA  HAT.  LV=L0WER  OBSERVED  VALUE. 

1  UV=UPPER  OBSERVED  VALU£,  OV=BINOMIAL  OBSERVED  VALUE 
CLS: PRINTP85."Do  you  wish  to  find  your  actual" 
PRINT0131,  AC.PHA  HAT  VALUE?" 

PRINTP202 , "Just  input  1 Y'  for  Yes  or  1 N 1  for  No" 

BS=INKEY$: I F  B$=""G0T0  13260 
IF  B$="N"  GOT 0 1 3100 
IF  B$="Y"  GOT014000 
GOTO  13240 

Y****  CALCULATION  OF  ACTUAL  ALPHA  HATS 
ON  QT  GOTO  14040,14300.14040,13100,14040 
CLS: PRINT@85,"Please  input  your  observed  1 T 1 " ; 

INPUT  B 

IF  B^O  OR  B>N  GOTO  14040 
IF  S>1  THEN  14620 
IF  T=1  THEN  14200 
IF  T=2  THEN  AT=CDF! (B+l) 

IF  T=3  THEN  AT=1-C0F! (B+l) 

CLS: PRINT086 , "ALPHA  HAT  is  equal  to";AT 
GOSUB  13020 
GOTO  13100 

LV=CDF! (B+l): UV=1-CDF!  (B+l) 

IF  LV<=UV  THEN  AT=2*LV  ELSE  AT=2*UV 
IF  AT>1  THEN  AT=1 


equal  to";AT 


=2*LV  ELSE  AT=2*UV 


GOTO  14180 

CLS: PRINTP83 ,  "Please  input  your  * TI '  value' 
INPUT  Tl 

CLS: PR I NT083 ,  "Please  input  your  * T2 '  value' 

INPUT  T2: CLS 

IF  T2>T1  THEN  14560 

IF  T=1  THEN  14480 

IF  T=2  THEN  AT =CDF! (Tl+1) 


•Jim  j  AAVl/tili 


-J  » J  *-11 


14460 

14480 

14500 

14520 

14540 

14560 

14565 

14580 

14600 

14620 

14640 

14660 

14670 

14680 

14700 

14720 

14740 

15000 

15020 

15040 

15060 

15080 

15100 

15120 

15140 

15160 

15180 

15200 

15240 

15260 

15270 

15280 

15300 

15320 

15340 

15360 

15380 

15400 

15420 

15425 

15440 

15460 

15480 

15500 

15505 

15510 

15515 

15520 

15540 

15560 

15580 

15600 

15605 

15620 

15640 

15660 

15680 

15700 

15720 

15740 

15760 

15780 

15800 

15320 

15840 

15860 

15380 

15900 

15920 


KseUt^ 


GOTO  14180 

LV=CDF! (Tl+1): UV=1-CDF!  (T2+1) 

IF  LV<=UV  THEN  AT=2*LV  ELSE  AT=2*UV 
IF  AT>1  THEN  AT=1 
GOTO  14180 

PRINTP83. "Please  check  your  numbers  and" 

PRINT@12&," re- input  ^TPand  'T2rfi 

GOSUB  130^0 

GOTO  14300 

IF  T=1  THEN  14680 

IF  T=2  THEN  AT=CDF!  (B+l) 

IF  T=3  THEN  AT=CDF!  (N-B+l) 

GOTO  14180 

LV=CDF! (B+l): UV=CDF! (N-B+l) 

IF  LV<=0V  THEN  AT=2*LV  ELSE  AT=2*UV 
IF  AT>1  THEN  AT=1 
GOTO  14180 

'  ****  INPUT  PARAMETER  SUBROUTINES 
'****  SUBROUTINE  FOR  N 

CLS: PRINT@125. "Please  input  N": 

INPUT  N: IF  N>6  THENGOTO  15120  ELSE  GOTO  15240 
CLS: PRINT0125  "Try  a  positive  number":  FOR  X=1  TO  1500:  NEXTX 
CLS:  GOTO  15040 
RETURN 

'****  SUBROUTINE  FOR  P 

CLS: PRINT0125, "Please  Input  P(robabi 1 ity)"; 

INPUT  P 

IF  P<0  OR  P>1.0  GOTO  15260 
RETURN 

CLS: PRINT@125,"P  must  be  between  0  and  1.0": FOR  X=1  TO  1500 
NEXTX 

CLS: GOTO  15160 
RETURN 

'****  SUBROUTINE  FOR  INPUTTING  ALPHA 

CLS:  PRINT@125, "Please  input  the  ALPHA  VALUE1'; 

INPUT  AT 

IF  AT<0  OR  AT>1.  0  GOTO  15420 
GOTO  15460 

CLS.  PRINTP125. "ALPHA  must  be  between  0  and  1.0" 

FOR  X=1  TO  1500: NEXTX 
CLS: GOTO  15320 
RETURN 

'*****  SUBROUTINE  FOR  CRITICAL  VALUE  TAIL  LOCATIONS  -T 
CLS: PRINT@83,"Please  select  your  critical  region" 
PRINTTAB(IO);  "1  -  Two  Tailed  Test" 

PRINTTABf 10); "2  -  Lower  Tail  Only" 

PRINTTAB( 10); "3  -  Upper  Tail  Only" 

PRINT@284, "PLEASE  MAKE  YOUR  SELECTION  <l-3>"; 

INPUT  T 

IF  T<1  OR  T>3  GOTO  15600 
GOTO  15620 

CLS: PRINT@126."T  must  be  a  1,2  or  3" 

FOR  X=1  TO  20£)0:  NEXTX:  GOTO  1$500 
CLS: RETURN 

r****  TH£  CALCULATION  OF  THE  BINOMIAL  PDF 

*********  IDENTIFICATION  OF  VARIABLES  ******.*** 

' P=PROBABI LITY 

' POF=PROBABI LITY  DENSITY  FUNCTION 
1 CDF=  CUMULATIVE  DISTRIBUTION  FUNCTION 
' N=THE  NUMBER  OF  OBSERVATIONS 
1 K=THE  NUMBER  OF  CLASS  ONE  OBSERVATIONS 
***********  BINOMIAL  PDF/CDF  CALCULATIONS  ********** 

CLS:  PRINT0128 , " I rm  calculating  CDF=0 
FOR  K=0  TO  N 
I=K+1 
PR0D=1 
M=K 

FOR  X=N  TO  (N-K+1)STEP-1 
IF  K=0  GOTO  16000 


********** 


*********** 


.*  ,  • .  • .  •  *  *v.»  * 


15940 
15960 
15980 
16C00 
16020 
16040 
16060 
16080 
16120 
16140 
16160 
16180 
16200 
16220 
16240 
16260 
16280 
16300 
16320 
16340 
16360 
16380 
16400 
16420 
16440 
164  '0 
164oJ 
16500 
16520 
16540 
16560 
16580 
16600 
16640 
16665 
16670 
16630 
16700 
16720 
16740 
16760 
16780 
16800 
16820 
16840 
16860 
16880 
16900 
16920 
16940 
16955 
16960 
16930 
17000 
17020 
17040 
17060 
17080 
17100 
17120 
17140 
17160 
17180 
17200 
17220 
17240 
17260 
17275 
17280 


PROD=PROD*(X/M) 

M=M-1 

NEXTX 

NEXTK 


RETURN 

i ****  < 

T  G( 


98. 


ROUTINE  FOR  ALPHA  VALUES  FOR  BINOMIALLY  DISTRIBUTION 
16280,16780,17140 

- ST 


TWO 


1 Hl-LOWER  ALPHA  OU 
’H2-UPPER  ALPHA  OU 
■HT-TOTAL  ALPHA  OU 
1 LN-"N"  -  LOWER  TA 
1 UN-"N"  -  UPPER  TA 
Al=(AT/2):A2=(AT/2) 
FOR  1=1  TO  (N+l) 


EC 
AE 

PUT  VALUE 
PUT  VALUE 
PUT  VALUE 
L 
L 


'LOWER  ALPHALOOP 


GOTO 
GOTO 
THEN 
•  UPPE 
GOTO 
GOTO 
THEN 


16460 
16520 
NEXT  I 
ALPHA 
16640 
16580 
NEXT  I 


F  CDF! 

J  CDF' 

IF  CDF! 

FOR  1=1 
IF  COF! 

IF  CDF! 

IF  ‘  ‘ 
rll=CDF! 

IF  LN< 

GOTO  16380 
H1=CDF!  (1-1):  LN 
■"  LN<0  THEf1 
16380 

1— CDF!  (1-1)): UN=I: HT=H1+H2 
T->4  GOTO  16670  ELSE  RETURN 
1 ~ CDF! (I)): UN= I  —  1 : HT=H1+H2 

T=4  OR  Qf= - 

B  18740 


LOOP 


IF 
GOT 
H2= 

IF 
H2= 

IF 

GOS _  . 

PRINT7AB(2)"1. 

PRINTTAB(2)"2. 

RETURN 
•  *  *  #  * 

IF  S=2 


.  :I-2 
LN=0 


2  THEN  16720 


1  T ' 
T' 


Is  less  than  or 
is  greater  than 


equal  to  LN 
S;UN:  PRINT 


LOWER  ONE  TAILED  TEST 
THEN  AT=AT/2 


N+l) 

-AT  GOTO  16860 
=AT  GOTO  16900 
) < AT  THEN  NEXTI 
)): LN=I-2 
'  LN=0 


1 


FOR  1=1  T 
IF  CDF! 

IF  CDF! 

IF  COF! 

HT=(CDF!  ( 

IF  LN<0  THEN 
GOTO  16955 
HT=( CDF' ,T-1)):  LN 
IF  LN--Q  HEN  LN=0 
IF  QT=4  OR  QT=  2  THEN  17100 
IF  5=2  THEN  HT=HT+HT 
IF  S=2  OR  S=4  THEN  UN=N-LN 
GOSUB  18740 

ON  S  GOTO  17060,17040,17060,17080 
PRINTTAB(4);  i1TT'  '  -  — +  — 

PRINHAB(4); ’"T1 

■"T1  is  greater  than  or  equal  to  ";(N-LN) 


PRINTTAB 
RETURN 

I  *  *  *  * 


is  greater  than  or  equal  to  ";(N-LN) 
is  less  than  or  equal  to  , LN: RETURN 

r  n  m  a  i  ^  H  An  i  i  ^  1  f  a  N  ,  (  K I  •  I  Kl\ 


ONE  TAILED  UPPER  TEST 
FOR  1=1  TO  (N+l) 

IF  CDF! ( I )>1-AT  GOTO  17220 
IF  CDF! ( I )= 1— AT  GOTO  17260 
IF  CDF! ( I )<1-AT  THEN  NEXT  I 
HT=( 1-CDF:  ( I ) ):  UN=I-1 
GOTO  17275 

HT=( 1-CDF! (1-1)):  UN=I 
IF  QT=4  OR  QT=2  THEN  17340 
GUSUB  18740 


S7 


7320 
7340 
8000 
3020 
3040 
3060 
8080 
3100 
3120 
18140 
18160 
13180 
18200 
18205 
13220 
18240 
13260 
13230 
18  300 
18320 
13340 
13360 
13330 
10400 
18420 
18440 
18460 
13480 
1 35  jO 
13520 
18540 
18560 
13580 
13600 
13620 
1064Q 
18660 
10630 
137'i0 


PR  I  NTT  AB( 5 ) ; T 
RETURN 

I  *  *  ★  ft 

GOSUB  15020 
GOSUB  15140 
GOSUB  1532" 
GOSUB  154| 
RETURN 

I  *  ft  *  * 


1  1 s  greater  than  UN 

SUBROUTINE  INPUTS 
'NUMBER 
'PROBABILITY 
'ALPHA 

'NUMBER  OF  CRITICAL  REGIONS  AND  LOCATIONS 


SUBROUTINE  FOR  DATA  INPUT 

Input  the  value  of  the  ";K,"TH  case":  PRINT 

INPUT  0:  DffK)=D:CLS:  PRINT0127 , "The  value  has  been" 

PRINT(3167,  assigned  to  case  '";K,"Tn 
PRINT:  PRINTTA8(7^("Is  this  correct  ' Y '  or  ' N ' ?" 


FOR  K=1  TO  N 
CLS:  PR  NT: PRINT" 
‘  INTTAB(  18): 


PR 


BS=INKEY$: IF  B$="n  GOTO  18240 
IF  BS-  •  " Y"  GOTO  18160 

SUBROUTINE  TO  DO  SHELL  SORT 


NEXTK 

I  *  *  *  A 


A=  1  442705042 
L=INT( LOG(N)*A) 
M=N 

FOR  11=1  TO  L 
M= INT(M/2) 

K=N-M 

FOR  12=1  TO  K 
1  =  12 
l1=I*M 


G0T0 ,861° 

PRINT 

NEXT  1 2 
NEXT  1 1 

1 FCR  X=1  TO  N:  PRINT  O' (X): NEXTX 
RETURN 

'  (LI) 


1  372U  G0T018540 

13740  r*  **  SUBROUTINE  PRINTING  OUTCOMES 

1  3 7 6 0  CLS:  ^RINTp85."  ALPHA  is  equal  to  " ; HT 

1 8  7 r  0  P  P  I N  T  T  AB(  S  )  '’reject  your  null  hypothesis  if":  RETURN 


SK 


APPENDIX  F 

CHI  SQUARE  DISTRIBUTION  PROGRAM 


1  CIS 


20  DIM  X»(15.15),C!(15),R!(15),N(3),N!(3) 
40  CLS:  PRlNT$123,"  Initial  program  load  in 
SO  A=8: N=0: CS=3:  t=l: MN=0:  SD=1:  TS=0 
60  FOR  1=1  TO  15: R! ( I )=0: NEXT 
80  FOR  J=1  TO  15: C! iJ j=0: NEXT 
100  FOR  1=1  TO  15: FOR  J=1  TO  15 


progress 


100  FOR  1=1  TO  15: FOR  J=1  TO  15 

110  X! ( I ,J)=0 

120  NEXT:  fsJEXT 

140  •***«*  VARIABLE  LIST 

160  '  X! (RT  CT)-ARRAY  OF  TABLE  CELLS//  R! (RT)-VECTOR  OF  ROW  TOTALS// 
162  '  C! ( CT } -VECTOR  OF  COLUMN  TOTALS//DF-QEGREES  OF  FREEDOM// 

164  '  EX-EXPECTED  CELL  FREQUENCY//  TS-TEST  STATISTIC// 

166  1  OT=  TEST  SELECTED//T=  TAIL  LOCATION  OF  CRITICAL  REGInN 
180  '  CI-CONFIDENCE  INTERVAL//  R .M-VARIABLES  IN  TOLERANCE  .iMITS// 
182  1  PA-THE  ALPHA  LEVEL//  Q-CERtA I NTY  LEVEL 


GREES  OF  FREEDOM// 
STATISTIC// 


REGION 

ANCE  wIMITS// 


166  1  OT=  TEST  SELECTED//T=  TAIL  LOCATION  OF  CRITICAL  REGInN 
180  '  CI-CONFIDENCE  INTERVAL//  R .M-VARIABLES  IN  TOLERANCE  .iMIT 
182  1  PA-THE  ALPHA  LEVEL//  0-CERTAINTY  LEVEL 
200  CLS: GOSUB  10000  ’THIS  IS  THE  CHI  SQUARE  MENU  SUBROUTINE 
220  CLS: PRINTP83 . "Do  you  wish  to  check  the  assumptions" 


UU.  V  W  I  (  V.-WV  «  WW  J  VU  f*  •  J  I  I  UW  v.  I  M.  t.  trUt  U  J  V  I  VM  J 

230  PRINTP123,  of  the  test  that  you  have  selected?" 

235  PRINT02O3 , " i nput  'Y'  for  yes  or  ‘N  for  no; 

240  E$=INKEY$: IF  B$=""  GOTO  240 
260  IF  B$-  Y"  OR  B$="y"  GOTO  320 
280  IF  BS="N"  OR  B$="n"  GOTO  900 
300  GOTO  240 

320  GOSUB  11000  'SUBROUTINE  OF  ASSUMPTIONS 
900  ON  QT  GOTO  1000,2000,3000,4000,5000,6000 

1000  '****  TOLERANCE  LIMITS 

1001  1 

1020  CLS: PR  I NT082 . "What  Confidence  Interval  do  you  want?" 
1030  PRINT:  PRINT(3l67,"(.  90,.  95,Your  Choice)"; 

1040  INPUT  CI:PA=1-Cl 

1060  CLS: PRINTP82, "What  is  the  population  proporti on?" 
1062  PRINT: PRINT@i70,"(Q  Value)"; 

1080  INPUT  Q 

1100  CLS:  PRINT(3127, "What  value  for  '  R ; 

1120  INPUT  R 

1140  PRINT0167, "What  value  for  ' M ; 

1160  INPUT  M: C ES : DF=2*(R+M) 

1170  IF  OF  =  0  THEN  1100 


1180  GOSUB  17000 

1200  N=.  25*TS*((l+Q)/n-Q)H.  5*(R+M-1) 

1220  PRINT086 . "To  be  h;CI;H  confident  you  need 
1230  PRINT012/ , "  an  N  >  or  =  to  ";N 
1240  PRINT: GOTO  5200 

2000  '****  MCNEMAR  TEST  FOR  N  >20 

2001  1 

2020  RT=2: CT=2:  DF=1 
2040  GOSUB  15000 


|+. 5*( R+M-l ) 
confident  s 


3SUB  15000 

2060  IF  (X! (1.2)+X! (2.1) )<20  THEN  2200 
2080  T'=( ( X! (I.2)-X! (2 , l))-2)/( X! ( 1 ,2)+X!  (2,1)) 

2100  Gl  %  160(50 
2120  GO,  5140 

2200  CLS;,  <INT@82, "Sorry  N  is  too  small  -we'll  start  over" 

2220  FOR  X=1  TO  1(500:  NEXTX:  CLS:  GOTO  40 

3000  '****  COCHRAN  TEST  FOR  RELATED  SAMPLES 

3001  ' 

3020  CT=0: RT=0: N=0: TS=0: CLS: PRINT0123 , "How  many  columns  do  you  have" 

3040  INPUT  C 

3060  FOR  1=1  TO  C 

3080  CLS: PRINTP123, "Input  the  column  "; I; "  total"; 


<20  THEN  2200 
2 )/( X! ( 1 ,2)+X!  (2,1)) 


3100  INPUT  K 

3120  CT=  CT+( K*K) 

3140  N=N+K 

3160  NEXT  I 

3180  TT=( C*(C-l)*CT)-( (C-l )*(N*N) ) 

3200  CLS:  PRINT&123 , "How  many  rows  do  you  have  ; 

3220  INPUT  R 

3240  FOR  1=  1  TO  R 

3260  CLS: PRINTP123, "Input  the  row  ";I;"  total"; 

3280  INPUT  K 

3300  RT=RT+(K*K) 

3320  NEXT  I 


( C*N)-RT 
>QB)):DF=(C-1) 


3340  CLS:  TB=(C*N)-RT 
3360  TS=(fTT)>(TB)):DF=( 
3380  GOSUB  16000 
3400  GOTO  5140 

4000  Y***‘  MEDIAN  TEST 

4001  ' 

4020  GOSUB  13500  ' INP 


4040  GOSUB 


INPUT  ROW  AND  COLUMN  TOTALS 
INPUT  TABLE  VALUES 


4060  I0=ABS( ( R!  m)-(R!(2))} 

4080  IF  ID<=1  THEN  4100  ELSE  4140 

4100  GOSUB  16000  'CALCULATE  CHI  SQUARE  DISTRIBUTION 
4120  GOTO  5140 

4140  CLS: PRINT043 , "Row  1  and  row  2  totals  should  be" 

4145  PRINTP83  "approximately  equal.  Do  you  want" 

4147  PRINT012S . "to  reinput  your  numbers?" 

4148  PRINT:  PRIATP203, '‘input  a  1 Y 1  for  yes  or  ' N '  for  no"; 

4160  INPUT  B$ 

4180  IF  B$="Y"  GOT04000 
4200  IF  B$="N"  GOT04100 
4220  GOTO  4140 

4999  Y 

5000  '***»  CHI  SQUARE  TEST 
5002  1 

5020  GOSUB  13500  'SUBROUTINE  FOR  COLUMN  AND  ROW  INPUTS 
5040  IF  (RT=2)  AND  (CT=2)  THEN  5500  'Gets  critical  tail  values 
5060  GOSUB  15000  ^SUBROUTINE  TO  FILL  TABLE 
5080  IF  CS=0  THEN  5700 

5100  GOSUB  16000  'SUBROUTINE  TO  CALCULATE  CHI  SQUARE  DISTRIBUTION 
5120  IF  T>1  THEN  5780 

5140  CS! =CS: T S J =TS: CLS: PRINTP47 , "Your  test  statistic  is  " 

5150  PRINT093, "equal  to  rt;TS! 

5160  PRINT0126 , "  your  YAlpfia  Hat'  i  s";  1-CS!  :  PRINT 
5200  GOSUB  13O0O 
5220  GOSUB  13100 

5500  '****  Selection  of  Critical  Region 

5520  CLS: PRINT@83 , "PI  ease  select  your  critical  region  " 

5525  PRINT013O,  1  -  Two  Tailed  Test" 

5530  PRI NTO 170, "2  -  Lower  Tail  Only" 

5535  PRINT0210 , "3  -  Upper  Tail  Only" 

5540  PRINT0284, "PLEASE  MAKE  YOUR  SELECTION  <l-3>"; 

5560  KS=INKEYS:IF  K$=tl"GOTO  5560 
5580  T=VAL(KS) 

5600  IF  C T< 1 j  OR  (T>3)  GOTO  5560 
5620  GOTO  5060 
5666  END 

5700  '****  Results  in  wrong  tail 

5720  CLS: PRINT083 , "Your  results  lie  in  the  wrong  tail." 

5725  PRINTTAB(8)"You  will  never  reject  your" 

5730  PRINTTABC I3)"nul 1  hypothesi s. " 

5740  GOSUB  13000 

5760  GOSUB  13100 

5761  Y 

5762  ' 

'★***  r  alrnla^irtn  f  a  i  1  a  1  r»Ka  wa  1  1  iq 


5780  '****  Calculation  of  one  tailed  alpha  valu 
5800  X9=TS:G0SU8  16140  Normal  calculations  qi 


5820  GOTO  5140 

cnnn  1****  n 


I  ue 

given  T>1 


GOODNESS  OF  FIT  TEST  -  CHI  SQUARE 


90 


i 


mww 


P 


[* 

f-r 

ry 

s 


6020 

6021 

6022 

6040 

6060 

6070 

6080 

6100 

6110 

6120 

6180 

6200 

6220 

6240 

6260 

6280 

6290 

6300 

6320 

6340 

6360 

6380 

6385 

6400 

6410 

6415 

6420 

6425 

6430 

6440 

6460 

6480 

6500 

6510 

6520 

7000 

7010 

7020 

7022 

7024 

7026 

7028 

7030 

7040 

7060 

7070 

7030 

7090 

7100 

7110 

7120 

7140 

7160 

7180 

7200 

7210 

7220 

7240 

7260 

7280 

7300 

7310 

7320 

7340 

7360 

7370 

7380 

7400 


classes"; 

UNIFORM,  EXPONENTIAL 
)f  yoi 

the" 


GOSUB  9000  „ 

CLS: PRINT082,  Please  input  the  number  of 
INPUT  CT 

ON  TO  GOTO  6060.7000,7000  'NORMAL,  _  , 

CLS: PRTnT@86. "Please  Input  the  mean  of  your 
PRlNT@126,  theoretical  distribution" 

INPUT  MN 

PRINt@166,"Now  enter  the  variance  of 
PRINT02O6," theoretical  di stribution" 

INPUT  VA 
FOR  J=1  TO 
PA=( 1/CT)  . 

GOSOB  18000 
N! ( 1 )=MN+SQR(yA}*IZ 
IF  J=1  THEN  o500 
CLS:  PRIN  ‘ 

PRINT012 
INPUT  C! 

N=N+C! (J 

N'^I 

"Jngutj  the  #  of  observed  data  points" 

CT): N^+Cf  Jit) 

@122, "How  many  parameters  did  you  estimate" 

1,_QR  2J"; 


IF J=1  THEN  6500 

:  PRINT@80 , "Input  the  #  of  observed  d 
NT@122 .  "are  >";  N!  (0);  “and  <  or  =fi;N! 


ata  points  that" 


CLS: NEXT 
PRINT082. 

PRINT@13< 

INPUT  C! { 

CLS:  PRINT . 

PRINT@215.  v„.  w. 

INPUT  XX:lF  X*<0  OR  XX*2  GOTO  6410 
DF=CT-(XX+1): EX=N/CT 
FOR  1=1  TO  CT 

TS=TS+( ( C!  (I ) -EX )  -*2 ) /EX :  NEXT 
GOSUB  16000 
GOTO  5140 

CLS: PRINT@83 , "Input  the  #  of  observed  data  points" 
PRINT0135  "<  or  -  ":N!(1); 

INPUT  C! ( 0 ): GOTO  63^0 

' ****  UNIFORM  DISTRIBUTION  AND  POISSON  DISTRIBUTION 
IF  TD=2  THEN  7100 
CLS: PRINT@45 ,  01 Q 


PRINTTAB 
PRINTTAB 
PRINTTAB 
PRINTTAB 
PRINTTAB 
E$=INKEY. 
IF  E$<>"M 


,  _  _  _ sses  are' 

.assumed  to_ run  from  0  to  ; ( CT- 1 ) 

i  An*' 


your  class* 
“  to"; ( 
assumption 


.  r ,  Please  note  - 
6)" 

61  "If  you  fai  1  this 
6)"to  return  to  the  input' menu 
l2r  STRIKE  'My" 

51"  ANY  OTHER  KEY  TO  CONTINUE" 
:  IF  E$="*GOTO  7040 
h  THEN  7080 


CLS: GOTO  40 
CLS: PRINTP125, "Please  input  LAMBDA"; 

INPUT  MN 
FOR  J=  1  TO  CT 

IF  TD=2  THEN  I=J  ELSE  I=J-1 

CLS: PRINT@46, "Input  the  number  in  column  ": PRINT@135 , "COLUMN-  "; I; 

INPUT  C!  ( J):  fsl=N  +C!  (J):  CLS 

NEXT:  IF  TD=2  THEN  6410 

'****  ENTERING  POISSON  LOOP 

EX=EXP(-MN)*N: EX=EX+. 5: EX=INT(EX) 

IF  EX=0  THEN  7500 
TS=TS+(fC! (l)-EX)-2)/EX 
1=0: PROD=l 
FOR  J=2  TO  CT 
1=1+1: PROD=PROD*I 

EX=N*^(MN|-Ij*EXP^-MN))/PR0D:  EX=INT(EX+.  5) 

TS=TS+( (C! (J)-EX)-2)/EX 
NEXTJ 

CLS: PRINT@122 . "Did  you  estimate  lambda  from  your  data" 

PRINT@210,  "  'Yy  for  yes,  or  ' N *  for  no" 

INPUT  ES 

IF  E$=" Y"  THEN  DF=CT-2: CLS: GOTO  6460 


7440 

7500 

7510 

7520 


the  expected  value  in  the  cell1 
need  to  collapse  cells." 


GOTO  7360 

CLS:  PRINT082, "Sorry , 

PRINT@125 i s  0,  you 
.  GOSUB  520(5 
7540  CLS: PRINT@83  "I  have  collapsed  all  cells  larger" 
7550  PRINT@130 ."tftan  or  equal  to";  I 
GOSUB  130(50 
GOTO  7360 

r****  THEORETICAL 


7560 

7580 

9000 

9020 

9022 

9024 

9026 

9040 

9060 

9080 

9100 

9120 

0000 

0020 

0022 

0024 

0026 

0040 

0060 

0080 

0100 

0120 

1000 

1020 

1040 

1060 

1080 

1090 

1100 

1120 

1130 

1140 

1160 

1180 

1200 

1220 


_  DISTRIBUTION 

CLS: PRINT041 , "Sel ect  your  theoretical 
PRINT@128.1-  Normal  Distribution" 

PR  I N T@ 1 68 , "2-  Uniform  Distribution" 

PR  I NT@208 , "3-  Poisson  Distribution" 

PRINT@284 ! ,JPLEASEiiMAKE_YOyR_SELECTION  <l-3>"; 


di stribution" 


K$=INKEY 
TD=yAL(K 
IF  (TD>3 
RETURN 

I  *  *  *  * 


IF  K$  =""G0T0  9060 
OR  ( TD< 1 )  GOTO  9060 


PRINTTABf 

PRINT"1-T 


SUBROUTINE  MENU  FOR  CHI-SQUARE  DISTRIBUTED  TESTS 


18); "MENU" 

....  olerance  Li 
PRINT  3-Cochran  Test  ,  . 
PRINT"5-Chi-Square  Test 
PRINT:  PRINT- PLEASE  M  _ 

'  lP  K$-  "GOTO  10060 

'or  (QT<1)  GOTO  10060 


K£=INKEY: 
•6 


-McNemar  Test" 
an  Test 

,  w  Goodness  of  fit" 

Selection  <i-6>"; 


SUBROUTINE  FOR  THE  TEST  ASSUMPTIONS 
ON  QT  GOTO  11040,11200,11280,11340,11420.11480 

PI  Q  ■  DDT  MTftQ£  ^Vcmi  cq!  orf  oH  Tr>1  ov'anro  I  irrnf  ClC 


CLS: PRINT@86,"You  selected  Tolerance  Limits":  GOSUB  12800 
GOSUB  13000 

CLS: PRINT043  "1-The  X 1 s  are  a  random  sample" 

PR  I NT@83 , " 2-drd i nal  level  data" 

PRINT@203,"T0  CONTINUE  -  STRIKE  ANY  KEY" 

P R I NT@243 ,  I f  you  fail  the  assumptions  then" 

PRINTP283,"to  return  to  the  menu  input  ‘M, 

E|=INKEY$;IF  E$='™  GOTO  11140 
IF  E$<>  M  GOTO  900 
CLS:  GOTO  160 

CLS:  PRINT@85,"You  have  selected  the  McNemar  test": GOSUB  12800 
GOSUB  13000 

1240  CLS:  PRINTP42 . " 1- 1 N1  equal  or  greater  than  20" 

1245  PRINTP82 . "2-the  pairs  (X,Y)  are  mutually1' 

"  PRINTP12S,  1 ndependent" 

PRINT0162 . "3-Nominal  level  data  in  two  categories": PRINT 
GOSUB  130(50 

CLS: PRINT@123 , " I f  you  failed  the  assumptions  then" 

PRINT@163 , "to  return  to  the  menu  input  '  M'  .GOTO  11140 
CLS: PRINT^84,"You  have  selected  the  Cochran  test": GOSUB  12800 
GOSUB  13000 

CLS:  PRINT@82,"l-Randomly  selected  treatment  blocks" 

P R I NT@ 122 , "2-Dichotomous  treatment  outcomes  (C 
CLS:  PRINTP83,"You  have  selected  the  Median  1 
GOSUB  13000 

1380  CLS:  PRINT042, "1-Independent  random  samples" 

1385  PRINT082 ,"2-At  least  ordinal  level  data" 

P R I NT@ 122 , "3-A1 1  populations  have  the  same 
GOTO  1110(5 

CLS: PRINTP43 , "You  have  selected  the  Chi  Square  test": GOSUB  12800 
GOSUB  13000 

1460  CLS: PRINT042  "1-Independent  random  samples" 

1470  PRINT082 ,"2-Each  observation  is  classified  into" 

1472  PRINTTABf 4 V'only  one  row  and  one  col umn": PRINT: GOTO  11100 
1480  CLS:  PRINTP42 ,  You  have  selected  the  Goodness  of  fit":G0SUE 
1500  GOSUB  13000 

1520  CLS: PRINT043 , " 1-An  independent  random  sample" 

1530  PRINTP83  "2-Data  level  appropriate  to" 

1533  PRINT@12S.the  distribution  selected": GOTO  11100 


1247 

1250 

1255 

1260 

1270 

1280 

1300 

1320 

1325 

1340 

1360 


1390 

1400 

1420 

1440 


:u  cauiiciii  y  1 

;nt  outcomes  (0,l)":GOTO  11100 
the  Median  tesirt: GOSUB  12800 


median" 


GOSUB  12800 


92 


12800 

12820 

12999 

13000 

13001 
13020 
13040 
13060 
13100 
13110 
13120 
13140 
13160 
13180 

13499 

13500 

13501 
13520 
13530 
13540 
13560 
13580 
13600 
13620 
13622 
13640 
13660 

14999 

15000 
15002 
15010 
15200 
15220 
15225 
15240 
15260 
15230 
15300 
15320 
15340 
15360 
15380 
15400 
15420 
15440 
15460 
154?0 
155C0 
15520 
15523 
15525 
15530 
15540 
15545 
15543 
15550 
15560 
15580 
15600 
15620 
15640 
15660 
15700 
15705 
15710 
15720 
15740 
15760 
15770 


PR  I NTTABf  6  j  this  test  assumes  you  meet  the 
PRINTTA8C 10 ) '' fo  1  lowing  assumptions":  PRINT:  R 


"fol lowing  assumptions":  PRINT:  RETURN 
SUBROUTINE  INKEYS 


PR  I  NTT  ABf  7 ) •  "TO  CONTINUE  -  STRIKE  ANY  KEY"; 

KS=INKEYS:  if7  K$="*  GOTO  13040 

RETURN 

CLS: PRINTP83 .  "Strike  ‘C1  to  do  another  test" 

PRINTP120,  If  you  are  finished  strike  any  other  key 

KS=INKEY$: IF  K§-  "GOTO  13120 

IF  KS-C"  GOTO  13180 

CLS: END 

ijLS:GOTO  40 

J****  SUBROUTINE  ROW  AND  COLUMN  INPUTS 

CLS: P R I NT083 , "PI ease  enter  the  number  of  rows  in" 
PRINTP127,  the  contingency  table"; 

INPUT  RT 

PRINT(3163."  Now  enter  the  number  of  columns"; 

INPUT  CT: If  (QT=4)  AN0  (RT<>2)  THEN  13620 


INPUT  CT-  IF  (QT=4)  AN0 
DF=(CT-iWRT-1):  RETURN 
CLS: PRINTP125  "You  are 

DD  TMTTftD/  11  \  n  i 


CLS: PRINTpl25  "You  are  permitted  only  ‘21  rows" 
PR  I NTTABf 1 1  )"i n  the  Median  Test* 

GOSUB  13000 
^LS: GOTO  13520 

’****  SUBROUTINE  TABLE  INPUTS 

i 

N=0 

FOR  I  =  1  TO  RT:  FOR  J=1T0  CT 

CLS: PRINTP46." Input  the  number  in  the  cell  in" 
PRINTPl 35 , "ROW-  I: PRTNT0175. "COLUMN-  fi;J 
PRINT0251 ."THE  NUMBER  ft; 

INPUT  X’ fl.J): C^S 
NEXT:  NEXT 
FOR  1=1  TO  RT 
R!  ( I)=0 
FCR  J=1  TO  CT 

R!  ( I  )=R!  (I)  +  X!  (I.J) 

NEXTJ 

N=N+R!  (I) 

NEXT  I 

FCR  J=1  TO  CT 
C! (J)=0 
FOR  1=1  TO  RT 

C!  (J)=C!  ( J  )  +  X!  (I.J) 

NEXTI: NEXTJ 

IP  QT=4  GOTO  15770 
IF  QT=2  GOTO  15660 
IF  T-*2  GOTO  157CO 
IF  RT<>2  AND  CT^>2  GOTO  15560 

RETURN 

FOR  I  =  1  TO  RT:  FOR  J=1T0  CT 
EX=R! (I)*C! f J ) / N : TS=TS+((X! ( I , J)-EX)-2)/EX 
IF  EX<5  THEN  CF=CF+ 1 
IF  EX-1  THEN  C0=C0+1 
NEXTJ: NEXTI 
RETURN 


IF  (T=2)  AND  (TS-0)  THEN  CS=0 
IF  (T=3)  AND  (TS<0)  THEN  CS=0 
RETURN 

Sl=(N-2)/(R! (1)*R! (2)) 


THEN  CS=0 
THEN  CS=0 


s 3=1  c*  ( ; );R'  cn  n 

.s:=s>C( 


sidi  cm  n 


MS1*S2  RFCRN 


,r.“ 

16400 

1 6  5  C  0 
■  r  s  >  - 
1  c  5  4  3 
16560 
iced 

1 6630 
16620 
1  c  6  4  o 
1  c  660 
16999 
17:90 

1  7  n  •  o 

•  '  o  c 

170C4 
1  7040 
1  7041 
1 7060 

1 7  0  60 
171^0 
17100 
17140 
17160 
17130 
17200 
17220 
17240 
17260 
17999 
1SC00 
13002 
18011 
1S020 
18040 
18060 
13030 
18140 
18150 


;***•  SDRCCINE  C  lA.lv.J’E  CC-SQuARE  DISTRIBUTION 

I  r  ’3-0  T^EN  16999 

CS'”-'R:nMM^:'tt  caMat’rg" 

•r  •'r.  4-  OC^C  ’ 6 "R C 

r..;.  N  IRMA  MI  5*  RIB,*  II  ON  APPROXIMATION  STEPS 
V^(I('S  QPJM  3J-1W  C*DF))  SQR(2):  V9=Y9*3*SQR(  DF) 

I*  . ...  -4  5  i^N  J > - 1  jO ,  j  Ir4?0 
C  >  3  -  4  3  'HEN  2S=I ■  oJTQ  Ib4J0 
V'M  M’  x --TS  ts  =  vj 

. . ormal  : adulations 

XM3  MABSl  (  TS-MN)  SO)  Rl=EXP( -rS*TS)/2)/2.  506628 

MI  (  >.  332b'*TS):  t .  5-RL*((.  937298M-.  120 1676 ) *TS+.  4361836 ) *TS 

; h: vn_t-en  z=.  d-M  ts=xl 

'this  ENDS  NORMAL  CALCULATIONS 
T  S  =  X  9  ■  C  S=2 
recrn 
dmc:  C9=i 

D9=C9*0l:  0l=DL-2: IF  DL"=2  THEN  16520 
P9=:NT((DF-i).  2):  N9=EXP( P9*LCG( TS )-TS/2)/D9 
DL-C":  E9=. 0000001:  F9= 1 E+07:  L9=l: M$=l: J9=l 
IF  DF-:*!NT(DF  2 )= 1  THEN  J9=SQR( 2/TS/3. 141592654) 

0L=Dl+2: M9=H9*TS/DL:  L9=L9+M9 
IF  M9  •  F9  THEN  CSC: RETURN 
IF  M9'E9  GOTO  16600 
l9=l9-M9:  CS=J9*N9*  L9:  RETURN 
CS  =  0: RETURN 

;*«**  INVERSE  CHI  SQUARE  DISTR.  BY  RECURSION 
CLS: PR  I NTOl 30 ."I'm  calculating" 

Mi&mA 

GOSUB  16000 

IP  (1-CS)<PA  THEN  17120 
NjiO^N^1):  GOTO  17040 

CLS-"pRINT@134  "Still  Going":  IF  FL=3  GOTO  17240 
IF  FL=5  GOtO  17260 
IF  FL=  6  THEN  17220 
A=A/4: GOTO  17040 
CLS: RETURN 

PRINT0170,  "  I'm  nearly  fi ni shed": GOTO  17180 
PRINT017O,  "  The  last  loop  to  do": GOTO  17180 

;****  INVERSE  NORMAL  DISTRIBUTION 

P8=PA: IF  PA>. 5  THEN  PA=  1-PA 
TL=SQR(LOG(l/(PA*PA))) 

RL=(. 0l0328*TL+. 802B53 )*TL  +2. 515517 
RL=RL/(((. 001308*TL+. 189269)*TL+1. 432788 ) *T L+ 1 ) 

IZ=RL-TL: IF  P8>. 5  THEN  IZ=-IZ 

RETURN 

END 
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